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ABSTRACT 


Context. The third data release by the Gaia mission of the European Space Agency (DR3) is the first release to provide the community 
with a large sample of observations for more than 150 thousand Solar System objects, including asteroids and natural planetary 
satellites. The release contains astrometry (over 23 million epochs) and photometry, along with average reflectance spectra of 60518 
asteroids and osculating elements. 

Aims. We present an overview of the procedures that have been implemented over several years of development and tests to process 
Solar System data at the level of accuracy that Gaia can reach. We illustrate the data properties and potential with some practical 
examples. 

Methods. In order to allow the users of Gaia DR3 to best exploit the data, we explain the assumptions and approaches followed in the 
implementation of the data processing pipeline for Solar System processing, and their effects in terms of data filtering, optimisation, 
and performances. We then test the data quality by analysing post-fit residuals to adjusted orbits, the capacity of detecting subtle 
dynamical effects (wobbling due to satellites or shape and Yarkovsky acceleration), and to reproduce known properties of asteroid 


photometry (phase curves and rotational light curves). 

Results. The DR3 astrometric accuracy is a clear improvement over the data published in DR2, which concerned a very limited 
sample of asteroids. The performance of the data reduction is met, and is illustrated by the capacity of detecting milliarcsecond-level 
wobbling of the asteroid photocentre that is due to satellite or shape effects and contributes to Yarkovsky effect measurements. 
Conclusions. The third data release can in terms of data completeness and accuracy be considered the first full-scale realisation of 
the Solar System survey by Gaia. 


Key words. Solar System: minor planets - orbits - astrometry - photometry; Gaia mission 


1. Introduction surveys, fostering new science achievements 

2007). With the accumulation 
The potential capacity of Gaia to provide an outstanding survey Of more accurate information about the mission performance, 
of Solar System objects (SSO) became clear already during the applications of the Gaia Solar System survey ere identified 
preparation studies (Hestroffer et al.]1999 [Mignard]2002). The in the improvement of our knowledge of dynamical properties 
expected sample of ~3x10° objects, including astrometric posi- Of asteroids through the ultra-accurate astrometry 
tions, photometry, and spectra, and the unprecedented accuracy, |2008}|Bancelin et al.|2012), leading in particular to the determi- 
were soon considered unique properties in the landscape of large Nation of asteroid masses (Mouret et al.!2007), the measurement 
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Fig. 1: Positions of the transits of sources published in Gaia DR3 in a full-sky Hammer-Aitoff projection in equatorial coordinates. 
The origin of RA and Dec is at the centre, north is up, and a increases from right to left. Darkness is related to the density of 
observations. Some density patterns related to the Gaia scanning law are clearly present. It is interesting to note the emergence 
of the highest stellar density regions as less dark areas close to the Galactic centre (around RA=270° and Dec=-30°), where the 


efficiency of unambiguous identifications decreases. 


the 


of the Yarkovsky effect (Delbo et al.|/2008 2015), 


discovery and characterisation of asteroid satellites (Pravec & 


Scheirich|2012)|Oszkiewicz et al./2013), the improvement of the 
dynamical models of satellite orbits (Arlot et al.)2012), a long- 
standing impact on ground-based observations of stellar occulta- 
tions (Tanga & Delbo}2007), and new tests of General Relativity 
(Hees et al./2018}|Hestroffer et al./2009). Epoch brightness mea- 
surements and low-resolution reflectance spectra were identified 
as an unprecedented source of knowledge about asteroid phys- 
ical properties: global shape properties and rotation parameters 
(Cellino & Dell’Oro|2012), compositions, and taxonomic clas- 


sification (Delbo et al.}2012). 
While Gaia DR2 provided a very limited and preliminary 


high-quality sample of astrometric and photometric data 


Collaboration et al.|2018), Gaia DR3 for the first time reaches 


the level of quality, variety, and volume that was expected for 
the Solar System (Fig. [I). The goal of this article is to illustrate 
the properties of the processing pipeline and the quality of the 
data that are obtained through the example of some significant 
applications. 


The implementation of the data processing for SSO was a 
long process that originated in preliminary studies that started at 
the end of the 1990s. These led to a functional analysis of the 
possible pipeline in 2006, several years before the launch of the 
satellite. Over time, the different processing modules have been 
developed, tested, qualified, and gradually entered into opera- 
tions with each data release. Solar System objects benefit from 
the improvement of all aspects of the Gaia data processing, with 
an increase in data quality from one release to the next. At the 
same time, the structure of the pipeline increases in complexity. 
It treats new features and produces a more complete data set at 
its output. 


While the fast daily processing that feeds asteroid alerts was 


illustrated elsewhere (Tanga et al.|2016 2021), we 


focus here on the procedures that have been implemented for 
Gaia DR3, aiming to exploit the whole accuracy of the data. We 
also intend to show how Gaia data, which are peculiar in many 
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aspects, should be used in practice. By doing so, we illustrate 
their potential for science with the example of some applications. 
Conversely, processing and validation of asteroid spectra are not 


discussed here as they are extensively presented by 
Taboration] (2022). — 

The article starts with a summary of the peculiarities of Gaia 
observations for the Solar System and of the general properties 
of the data present in Gaia DR3 (Section [2). We then illustrate 
the data-processing pipeline (Sect. [3). starting with a descrip- 
tion of the input data. We then describe the principles adopted to 
identify SSOs in the general data stream, to derive the astrome- 
try, and to compute the orbits and the calibrated photometry. We 
provide information about the result of a match with a recent or- 
bit catalogue for moving sources that are listed as not identified 
in Gaia DR3 (Sect./4). The quality of the astrometric and photo- 
metric data is then illustrated by several examples in Sect. [5]and 


(6) 


2. Solar System data in Gaia DR3 
2.1. Summary of the general properties of Gaia observations 


We recall some basic properties that strongly drive the process- 
ing, the results, and the use of the data with the relevant ter- 
minology. Exhaustive descriptions can be found in the 
documentation and in|Gaia Collaboration et al,|(2016). General 
properties of Solar System observations were also provided with 
Gaia DR2 in (Gaia Collaboration et al-[2018). 

The continuous rotation of the satellite results in the drift 
of all sources across the CCD matrix on the focal plane. Each 
focal plane passage (called transit”) can provide nine positions 
at most (in the Astrometric Field instrument, AF) and two low- 
dispersion spectra (in the Blue and Red Photometers, BP and 
RP). The AF is unfiltered and produces the G-band photometry. 
At the beginning of each transit, sources are first detected by 
the Sky Mapper instrument (SM). While the SM is essential for 
the on-board assignment of pixel windows to track the transiting 
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Fig. 2: Distribution of the observations of Solar System objects in solar elongation as a function of time. Colour is related to data 
density. Two overlapping sinusoids appear, corresponding to the variation in the directions in which the scanning plane of Gaia 
intersects the ecliptic due to the precession of the spin axis of the satellite. The peaks correspond to periods in which the ecliptic 
is crossed by the scan both at ~45° and in the opposite direction ~135°. In these conditions, the scanning plane is perpendicular 
to the ecliptic. Scattered data correspond to detections of objects at high ecliptic inclination. The nodes of the sinusoids are around 
quadrature (90° elongation), when the scanning plane cuts the ecliptic at 45°. An initial period without precession (enforcing the 
ecliptic pole scanning law) is visible. Three gaps appear, corresponding to technical operations preventing the normal collection of 
data. In the right panel, the cumulative distribution shows the preferential accumulation of observations at the extremes. 


sources, the astrometry that it provides has lower quality and is 
not published. 

The scanning motion of the Gaia telescopes combines the 
rotation of the satellite (period of six hours), the precession cy- 
cle of the rotation axis (68 days) on a Sun-centred cone, and 
the revolution around the Sun (one year). These three motions 
determine the typical timescales and locations when and where 
SSOs are observed. Single CCD measurements providing posi- 
tions during a transit are spaced by 4.4s. G-band photometry, 
averaged over the AF, combines observations over ~40 s. 

In Gaia DR3, the astrometry of an SSO is provided for each 
CCD. Hereafter, ’position” refers to this single CCD measure- 
ment. G—band photometry is also derived from CCD-level mea- 
surements, but is provided as an average value over all CCDs that 
are available during a transit. This is implicitly assumed when 
the brightness of an object is mentioned. In summary, astromet- 
ric data of SSOs are provided at individual CCD levels, while 
photometric data are provided as averages over a transit (see Ta- 
ble[Ip. 

The orientation of the scanning motion with respect to the 
Sun results in two avoidance cones with a semi-aperture of 45°, 
one centred on the Sun, and the other at solar opposition. As a 
consequence, observations are always obtained at solar elonga- 
tions in the range between 45°and 135°. The distribution within 
this interval is not uniform and strongly favours the extreme val- 
ues of elongation, where more time is spent by the scanning mo- 
tion (Fig. 2). 

The two telescopes on board Gaia sweep almost the same 
sky area 106 minutes apart (over a single rotation of the satel- 
lite). Conventionally, they are referred to as preceding and fol- 
lowing fields of view (FOVs), or FOV1 and FOV2. For a given 
asteroid, this scanning law usually results in short sequences of 
consecutive observations, separated by an absence of detection 
over several weeks or months. 

By design, focal plane pixels are rectangular with an aspect 
ratio of three. Their short side is oriented in the scan direction 
(along scan, AL). With the exception of the brightest sources 
(G<13), the signal is also binned in the across-scan (AC) direc- 
tion. Only the AL accuracy (at milliarcsecond level) is fully pre- 
served, while only an approximate position is available for AC 
(with an accuracy of about one arcsecond). To reduce telemetry 
volumes, only a window (a limited surface of pixels around each 
source) is acquired by the image-processing system on board 
Gaia. As the window coordinates are computed at the beginning 


of the transit following the detection in SM, and then propagated 
across the focal plane by the predicted motion expected for a 
stellar source, moving objects (e.g. the SSOs) will drift with re- 
spect to the window centre, and their signal can be truncated at 
the edge of the window. 

The global astrometric solution of Gaia, provides absolute 
positions on the sky. This also applies to SSOs. The exact defini- 
tion of the reference system and timescale is critical for the best 
exploitation of Gaia data (Sect.|2.2). 


2.2. Definition of the reference frames 


Because the astrometric positions resulting from the Gaia data 
are highly accurate, both the reference frame and the timescale 
must be clearly defined and well understood by the user. The in- 
dividual observations at the CCD level and the orbital elements, 
or equivalently, the state vector of the observed objects, must be 
distinguished first. 


2.2.1. Astrometric positions 


The positions given as Gaia-centric right ascension (RA) and 
declination (Dec) are derived from the local coordinates in the 
physical pixels and are transformed to astronomical coordinates 
as explained in Sect. |3.4.1} The final positions are given in the 
Barycentric Celestial Reference System (BCRS) with the origin 
at Gaia, and this is achieved by ultimately referring the attitude 
of Gaia to the Gaia-CRF3 axes, which are aligned to ICRF3 
(Gaia Collaboration, Klioner et al |2022). The directions pro- 
vided in the form of RA, Dec are similar to the astrometric posi- 
tions, meaning that they are corrected for the annual aberration, 
but not for the relativistic light deflection by the Solar System 
gravity field. This last correction for sources at a finite distance 
requires knowledge of the distance, which is known for all the 
SSOs in Gaia DR3, with the exception of a fraction of asteroids 
in the unmatched category. Taking the presence of light deflec- 
tion into account is left to users computing orbits, as for the light 
travel time. For the relativistic framework adopted for Gaia, we 


refer to|Klioner]|(2003} [2004). 


2.2.2. Orbital elements 


Orbits are computed from a least-squares fitting of a dynamical 
model to the observed Gaia-centric directions at transit times. 
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Fig. 3: General structure of the data processing. The turquoise (filled) boxes are modules developed specifically for Solar System 
objects, which run on the servers of one of the Data Processing Centers of Gaia (CNES in our case). White boxes represent input and 
output data. Yellow boxes are tasks that are developed, managed, and run by the Solar System scientific team that are essential to the 
preparation of the input data or to the validation of the output. Other modules providing input data, developed by other teams, are 
initial data treatment (IDT), intermediate data updating (IDU), the astrometric global iterative solution (AGIS), and the photometric 


pipeline (PhotPipe). 


This is detailed in Sect. [3.5] For each SSO, the unknowns are the 
six components of the state vector (position vector, velocity vec- 
tor) at a certain epoch, nominally, the median of the transit times 
of the SSO, in order to minimize the correlations. The state vec- 
tor is heliocentric and given in the ICRF axes. The Keplerian or- 
bital elements are derived from the state vector with the ecliptic 
as reference plane. In this context, the ecliptic is defined by two 
rotations from the ICRF, the obliquity and the origin bias. This 
results in a rotation matrix. The values used throughout the Gaia 
data processing are those from |Chapront et al.|(2002), with the 
origin bias of @piag = —55.42 mas defined as the right ascension 
of yicrs on the ICRF fundamental plane. In this definition, yjcrs 
is the intersection of the dynamical ecliptic with the ICRS funda- 
mental plane, which is not the same as the J2000 celestial equa- 
tor. The obliquity that is used corresponding to the inclination 
between the two planes at yicrs iS Gcrs = 84381'74110. These 
values are not exactly identical to those adopted in the SOFA 
standards (—52.928 mas and 843817412819). The yjcrs so de- 
fined is used as the origin of longitude in the ecliptic plane for the 
longitude of node in the orbital elements. It differs from yy2999, 
the intersection of the mean ecliptic with the celestial equator at 
J2000 by an angle of ~ 42 mas, with the sign convention such 
that this angle is the longitude of yj2999 referred to yjcrs. The 
Gaia-derived longitudes of node may therefore show a system- 
atic when compared to other sources of orbital elements using 
the origin of longitude in the celestial equator of J2000. 

Finally, the rotation matrix to transform a vector given in 
the ICRF to the same vector expressed in the ecliptic as defined 
above is 


R= RU, €icrs) RG, Ppias)- 


This matrix has been applied to the orbit determination to ex- 
press the state vector in the ecliptic frame before the heliocentric 
orbital elements were computed. The same matrix is used for 
the transformation of the covariance matrix, but in this case, the 
exact definitions are less critical. 
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2.2.3. Timescale 


Very early in the data processing, it was decided to use the TCB 
as the astronomical timescale for all the Gaia computations. This 
was a logical choice after the BCRS and the associated relativis- 
tic metric were chosen as the framework for the astrometric mod- 
elling. The on-board time tagging is calibrated against the TCB 
on the ground to obtain the correspondence between the two 
scales and to provide a final timing of all observations and Gaia 
events in TCB. Because of the unique accuracy of Gaia astrom- 
etry, the internal consistency was an essential requirement for 
the processing. This means that every ephemeris used in the pro- 
cessing has TCB as an independent variable for the Solar system 
(major and minor planets, natural satellites), but also for the or- 
bit of the spacecraft itself. This contrasts with the more common 
use of TDB for the public ephemeris and also as the timescale 
for the epoch of the orbital elements in astorb 
or at the Minor Planet Center. The transformation between 


the two scales is given by|Berthier et al.](2021) and|Klioner et al.| 


(2010) following IAU resolution 2006 B3'| 


TDB = TCB 
— Lg(JDrcpg — 2443 144.500 3725) x 86400 s 


~§55 10° 


where the time is expressed in seconds, and Lg = 
1.550519 768 x 10-° is a defining constant in the astronomical 
system of units. During the period covered by the Gaia DR3, the 
difference TDB — TCB is ~ —19s. 

As a purely indicative approximation, the UTC at the posi- 
tion of Gaia as derived from TCB is provided in the astrometry 
table gaiadr3.sso_observation of Gaia DR3. Nevertheless, 
UTC should not be used for an accurate exploitation of Gaia 
astrometry. 


' https: //iau.org/static/resolutions/IAU2006_Resol3. 
pdf 
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2.3. Population of Solar System objects 


In comparison to DR2, Gaia DR3 is more than ten times richer 
in terms of objects, a factor of >1.5 longer in time span, and pro- 
vides a more complete data set. A summary of the published data 
is provided in alana All Some general properties are illus- 
trated by Figs. [4]and All published data concern small SSOs, 
with the following main categories (Table|2}: main-belt asteroids 
(MBAs), near-Earth objects (NEOs), outer Solar System popu- 
lations, unmatched (unidentified) moving objects, and planetary 
satellites. 


Table 1: SSO data summary in DR3 at CCD level, at 
transit level, or per object. See text for explanations. 


Data type CCD-level transits 
astrometry 23 336467 3214776 
G-band photometry 3 069 170 


number of objects 


reflectance spectra 60 518 
orbits 154 741 
Table 2: Object types in DR3. 

Object type number of objects 
Atira 1 
Aten 43 
Apollo 230 
Amor 173 
Mars Crossers 1550 
Inner Main Belt 3305 
Main Belt 144 975 
Outer Main Belt 4940 
Jupiter Trojans 1550 
Centaurs 8 
TNOs 24 
Others 2 
Total asteroids 156 801 
Unmatched moving objects 1 320 
Planetary satellites 31 
Total 158 152 


3. Data processing for the Solar System 


All data products are the result of the data treatment pipeline de- 
veloped by scientists of the Gaia Data Processing and Analysis 
Consortium (DPAC), implemented in the computing facilities of 


the French space agency (CNES). Strict qualification and valida- 
tion protocols have been followed to ensure correct results and 
full consistency with the other subsystems of DPAC. 

The general structure of the data processing is illustrated in 
Fig. |3} We describe it below by focusing on the aspects that 
more strongly define the properties of SSO data found in Gaia 
DR3 to allow potential users in the scientific community to be- 
come familiar with the data and their exploitation. Orbits have 
not been computed by the core pipeline, but by an offline proce- 
dure (Sect.[3.5) with a dedicated validation. 

The core of the input data comes from the intermediate data 
updating (IDU) and from the astrometric global iterative solution 
(AGIS), wich are two components of the general data processing 
of Gaia. Details about their implementation principles can be 


found in (2016), {Rowell et al.](2020),|Torra et al.| 
(2020) and|Lindegren et al.|(2021). 


IDU in particular provides all the information required to re- 
construct the position and brightness of a source in its window, 
starting from the image parameter determination (IPD; i.e. the 
determination of the centroid of the signal). The data exploited 
for each transit are 


— the value of the magnitude determined on board by the video- 
processing unit (VPU), as a preliminary estimation. 

— The window class generated by the VPU. It is possible to 
reconstruct the window geometry for each strip based on this 
parameter. 

— The along-scan and across-scan window coordinates are re- 
lated to the timing and position of the window in the corre- 
sponding CCD (reconstructed by IDU). 

— The data also include the list of along-scan centroids x, 
(DU), 

— the list of the fluxes f, (photo-electrons per second, by IDU), 

— the list of the across-scan centroids y, for bi-dimensional 
windows (IDU), 

— and a list of flags, generated by IDU, which describes the 
quality of the IPD output and the encountered issues, such as 
the formal errors, the goodness-of-fit (GoF) to the PSF/LSF 
model, the presence of secondary peaks in the window, or 
the background estimation. 


Along with (x;, ys), the window reference system coordi- 
nates (WRS), the corresponding epochs are present, expressed 
in the internal on-board mission time line (OBMT; for its defini- 
tion, see|Gaia Collaboration et al.]2016). AGIS provides all the 
information required to calibrate the astrometry and to transform 
the WRS coordinates into BCRS. A particularly relevant input is 
a pre-computed list of transits that provides a first identification 
of observations associated with SSOs. We illustrate this in the 
following section. 


3.1. Selection of sources 


Two main reasons prevent an automated selection of the SSOs 
within their dedicated pipeline: first, the huge volume of data 
in the automated pipeline, of which only a small fraction must 
be selected and exploited; and second, the need to optimise the 
extraction of SSOs while avoiding contaminants as much as pos- 
sible. 

The SSO selection must therefore be performed by a pre- 
processor and is then exploited at the ingestion of the pipeline to 
select only the corresponding IDU and AGIS data. We adopted 
two approaches. The first approach is the most relevant for pro- 
cessing the bulk of the SSO sample. It is devoted to objects that 
can be identified by a direct match to their predicted positions. 
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The second approach retrieves sequences (bundles) of moving- 
object detections that are not matched to known SSOs. 

For the first selection, several criteria were applied to obtain 
a subset of objects that was about ten times larger than in the 
previous Gaia release (DR2). 


— We aimed for between 100000 and 150000 objects in the 
input list that were representatives of all the broad categories 
of asteroids were sought, such as NEOs, MBAs, Jupiter Tro- 
jans, and transneptunian objects (TNOs). Transits of several 
planetary satellites were also included. Comets alone are not 
present in Gaia DR3. 

— A transit was not selected if a star, another SSO, or a con- 
taminant generated by a bright star was found too close to 
the object during its observation by Gaia. 

— Each selected SSO had to be detected on at least eight transits 
over the 34 months covered by the Gaia DR3 data. 


Known SSOs were searched for by matching all the observed 
transits (from the output of the initial data treatment, IDT; 
to computed transits of SSOs over the Gaia 
DR3 time span. The computed positions as seen by Gaia were 
obtained by the available information on the position of Gaia in 
space (the satellite orbit), the scanning law, and a numerical in- 
tegration of the SSO motion. This last procedure starts from the 
osculating elements and osculating epoch given in the astort’| 
database (Bowell et al.]1993). As the selection had to be finalised 
well before the operations of the SSO pipeline, the astorb ver- 
sion we used is that of 13 December 2017. 

The SSO cross—matching proceeds in two steps: first by the 
crossing time (required match within 0.1 s), and then by the sky 
coordinates within a window of 1’’5. After the list was filtered 
for possible contaminants, the final input selection had 3 513 248 
transits for 156 837 known asteroids. To meet the criteria men- 
tioned above for Gaia DR3, a search for numbered asteroids was 
sufficient. Unnumbered asteroids do not appear in this release, 
with the exception mentioned further below. 

The search for observations of planetary satellites observ- 
able with Gaia has proceeded like for the asteroids. Their pas- 
sages were first predicted using the ephemeris of the satellites 
provided by V. Lainey from IMCCE. Due to the angular prox- 
imity of the planet, the number of contaminants is much higher 
than average, and the selection of a single match may become 
somewhat arbitrary at times. This is the case in particular for the 
Galilean satellites. Unlike for the asteroids, there has been no fil- 
tering for the number of detected transits. Out of the 44 possible 
satellites in the appropriate magnitude range (2 for Mars, 18 for 
Jupiter, 14 for Saturn, 7 for Uranus, and 3 for Neptune), 42 have 
been matched to at least one IDT transit. The two missing ob- 
jects are the Jovian satellites Amalthea and Thebe, whose signal 
is dominated by the noise due to their proximity to Jupiter. 

The astrometry of the potentially unknown moving sources 
constitutes a smaller data set that appears in Gaia DR3 for the 
first time. As mentioned above, their selection is based on an 
independent search that heavily relies on the along-scan motion 
(Vaz) that is reliably provided by IDT (see Sect. 5.2 in 
only after 1 December 2016. Therefore, for Gaia 
DR3, there were altogether only six months of data that could be 
exploited for this search. 

The selection recovers the position of sources above a chosen 
AL motion threshold (> 1.5 mas s~!) provided they can be paired 
with at least one other after an interval of time corresponding to 


* Ihttps://asteroid.lowell.edu/main/astorb/ 
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transits in the preceding-following field of view (or following- 
preceding), consecutive or separated by a satellite rotation (6 
hours; Sect. 2.1). The pairing is validated only if the estimated 
Vaz based on the two positions and the two Vay, values of the po- 
tential pair are also compatible. In a second step, all the pairs are 
examined in succession to connect to other pairs that may belong 
to the same source and are thus appended to form longer chains, 
with three or more observations in a sequence found in consec- 
utive IDT runs. In chains of four or more transits, the number of 
constraints is so large that the link of the detections to a single 
source is certain. 

The acceptance thresholds and, at the end of the process, the 
reliability of the detection, depend on the number of consecutive 
observations. While sources with only two or three transits dom- 
inate the sample, approximately 10% of the sources have four 
or more transits, the largest number of transits is 42. In the pro- 
cess, all the known asteroids (whose ephemerides are computed 
based on the astorb file mentioned above) are also found and 
discarded. 

The final input list of unidentified SSO for DR3 comprises 
4522 transits of unmatched asteroids, corresponding to 1531 
groups of chained transits. Unmatched sources do not necessar- 
ily correspond to new asteroids. This sample includes asteroids 
that were not in the orbital data at the end of 2017 because their 
orbit was not available at that time or because it was too poor to 
have a successful position match in the window of 1’’5. This is 
a vivid illustration that the population of asteroids is not sharply 
divided between known and unknown SSOs. These are just the 
two boundaries of a continuous spectrum of knowledge. A fur- 
ther exploration of this selected sample of unmatched objects, 
based on an updated version of the orbit data base, is provided 
in Sect.[4] 


3.2. Identification of Solar System objects 


The first step of the SSO pipeline is identifying each source 
entering the processing based on the average transit position 
provided by the initial data treatment (Sect. 6.4 in 
et al.}2016). This is basically done by matching the approximate 
source position (~0‘’1) to the ephemerides of SSOs. This proce- 
dure essentially duplicates the one adopted to build the input list, 
but with the fundamental difference that most of the potential 
sources that are not SSOs have not been ingested in the pipeline. 
The new identification is thus performed on a rather clean data 
set, in which sources that are not SSOs, in overwhelming number 
at the input list selection, are now reduced to a minimum. An- 
other relevant difference is as follows: at this level, no attempt is 
made to identify unmatched sources or to rebuild their bundles. 
As time passes, the inventory of known SSOs grows. Since 
the time of the Gaia launch, hundreds of thousands of new as- 
teroids have been discovered or their orbits were considerably 
improved, bringing the total number of known asteroids with 
an orbit to more than 1 170000 as early as 2022. A database 
of precomputed ephemerides of all known SSOs is regularly 
updated for the need of the processing of Gaia data. Time- 
dependent Gaia-centric positions of SSOs are arranged by using 
a HEALPix spatial index (Gérski et al.|2005), a grid resolution 
Neide = 2!°, and a time resolution adapted to_each object. They 
are stored into an Apache Cassandra databasq>| During each pro- 
cessing cycle, packets of transits are cross-matched with known 
SSOs extracted from the database. For this purpose, the pairs 
{HEALPix, Epoch} of each transit are used to extract a sample 


3 https://cassandra.apache.org/ 
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Fig. 4: Visual portfolio of some data properties for asteroids in Gaia DR3. In the top panel, we show the distribution on the semi- 
major axis of objects belonging to the main populations (eight Centaurs are excluded). The bin size is 1/30 au in the left part and 
1/ au for TNOs. Middle panel: Distribution of main-belt asteroids in the proper elements as provided by Astdys on the a, e (left) and 
a, sin(1) planes (right). Colours represent sin(i) and e, respectively. The left bottom panel presents the distribution of the published G 
magnitude per transit. The cut imposed at G = 21 is visible. The right bottom panel shows that the distribution of the H magnitudes 
of the objects, as provided by the Minor Planet Center, is strongly dominated by main-belt asteroids, with a contribution from NEOs 


for the faintest sources. 


of zero to a few dozen SSO candidates. Their Gaia-centric accu- 
rate positions at the transit epoch are then recomputed by means 
of two-body numerical integration perturbed by n bodies, pro- 
viding equatorial coordinates that can be directly compared to 
measured transit coordinates. 


The first criterion of candidate selection relies on the accu- 
racy of SSO orbits. For each target, the ephemeris uncertainty 
at the epoch of each transit is computed based on the lo RMS 


(o>) of its orbit, again adopting the astorb database. A candi- 
date is retained if its current ephemeris uncertainty (CEU) fulfils 
the condition 


CEU = 0, + (t-t) Go < €, (1) 


where ¢ is the observation epoch of transit, fo is the reference 
epoch of the orbital elements, c-, is the rate of change of o,, and 
€ is a given threshold. The adopted value e« = 10” leads to the 


Article number, page 7 of 34 


A&A proofs: manuscript no. Tanga_DPCP150_final 


10! 


data density 


10° 


Fig. 5: Distribution of H magnitudes for asteroids in the semi- 
major axis range from NEOs to Trojans. The trend followed by 
the detection limit is visible. In the mid-belt at 2.8 au, the limit 
G~21 corresponds to H~16.5. 


rejection of all SSOs with uncertain orbits, which could lead to 
spurious identifications. 

The second criterion takes into account the relative positions 
of SSO candidates compared to the recorded transit position, as 
described by [Pineau et al.](2011) in their cross-correlation algo- 
rithm. The SSO position 1s projected onto a 2D plane centred on 
the transit position, so that the relative coordinates of the SSO 
are x = d, y = O, where d is the angular distance between the 
two sources calculated by the haversine function, 


24 fOg 0 .9{@s—-a 

d = 2arcsin ,/sin? (=) + sin’ (5) cos 6,COS6s, (2) 
where @;, 6; and as, 6; are the equatorial coordinates of the tran- 
sit and the SSO candidate, respectively. In this plane, an SSO 
candidate is retained if its coordinates satisfy the condition 


g i (3) 


Ox, V1 — (co x,7y,)? 


where k = 3.43935 is the 2D completeness value for a 30 crite- 


rion (e.g. 99.7%), and where oy, = ,/o% +02 and ppox,y, = 


P10x,0y, + PsOx,7y, represent the uncertainties on the positions 
of the transit (t) and the SSO candidate (s) expressed by their 
covariance matrix, assuming Gaussian uncertainties (see Ap- 
pendix A of (Pineau et al.J2011). For SSOs, the positional uncer- 
tainty is taken as the current ephemeris uncertainty, for instance, 
Oo, = 05, = CEU. For transits, the positional uncertainty is fixed 
to 0.5’’, providing a large margin over the formal uncertainty of 
IDT (about 006 in each coordinate). 

A third criterion based on the difference in magnitude be- 
tween the observed transit and SSO candidates might be used 
in principle to distinguish between different candidates. Never- 
theless, the uncertainty on the predicted apparent magnitudes of 
many SSOs can reach values of about 1 mag or more because 
their albedos, light-scattering properties, and shapes are poorly 
known. No magnitude-based criterion was used in Gaia DR3. A 
consequence of this choice is that a faint source very close to an 
asteroid (e.g., a possible satellite) might in principle be matched 
to the asteroid itself. Despite the fact that in Gaia DR3 close 
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couples are filtered at the input list generation, we cannot totally 
exclude that such double detections exist. 

While planetary satellites have an identification based on the 
computation of their own ephemeris, no specific procedure is im- 
plemented to identify satellites of asteroids. As they share very 
similar coordinates, the satellites and the main body of the sys- 
tem can be given the same identifier. For instance, in Gaia DR3, 
the dwarf planet (134340) Pluto and its main satellite Charon are 
both identified as “(134340) Pluto”. 

If more than one SSO candidate satisfies the identification 
criteria, there is no obvious method to identify the correct object 
at this step of the processing. The best candidate is thus selected 
by calculating the quadratic distances between each observed 
transit and the corresponding SSO candidates. When there is 
more than one possible choice, the object minimising the dis- 
tance above is selected. With the grid resolution chosen, this sce- 
nario is fortunately very unlikely. 

The validation of the identification process has shown that 
the rate of correct identifications is very close to 100%, with 
an uncertainty smaller than 1%. This mainly comes from un- 
certainties on the positions of some SSO candidates and from 
the presence of some unfiltered contaminants. The identification 
process successfully recovers 99.97% of the transits in the input 
list. The small fraction of non—identified transits, negligible in 
practical terms, is mostly due to minor differences in the adopted 
ephemerides that come from differences between osculating ele- 
ments that are used to build the input list and those that are used 
to compute processing ephemerides. 

With the identification, a computation of the SSO 
ephemerides is performed for all the known sources. Some ancil- 
lary data, such as the distance and the apparent motion on the sky 
are computed, stored in an appropriate table, and are propagated 
to the pipeline where they remain available to other processing 
modules. 


3.3. Raw centroid processing 


The distribution of the collected photo-electrons inside a pixel 
window contains the raw information about the location of the 
source. The determination of the average position of any source 
from this distribution, the “centroid”, is carried out by IDU by 
means of a fit to a suitable model (Rowell et al./2020). The model 
is based on two assumptions: (1) the source is point-like, and (2) 
its image on the focal plane moves in the along-scan direction at 
a rate that exactly matches the charge transfer in the CCD. The 
free parameters of the model are the mean position of the source 
(centroid) and its intensity (flux). 

In the case of SSOs, hypothesis (1) is fulfilled in the very 
large majority of observing circumstances because only larger 
objects (order of thousands) can cause detectable signal distor- 
tions with respect to the point spread function (PSF) of the Gaia 
instrument. Conversely, assumption (2) is in general not valid 
for SSOs because the proper motion of SSOs with respect to 
stars produces a systematic shift of the photo-electron distribu- 
tion with respect to the scan motion rate. 

While a pure shift without smearing entails only a different 
but correct value of the centroid, the smearing introduces sys- 
tematic biases both in centroid and flux determination because 
the signal is truncated at the window edge. The magnitude of the 
bias increases with the distance of the centroid from the centre 
of the window and depends on the velocity of the source. Details 
are given in the|online documentation) 

In order to mitigate the impact of the centroiding bias on the 
final data quality, IDU positions are rejected that due to their 
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proximity to the window edges are expected to have a bias ex- 
ceeding the formal error on centroiding. This filtering is per- 
formed just before the astrometric processing module discussed 
in the next section. As the effect of the shift accumulates along 
the transit, the last AF columns have a stronger rejection prob- 
ability. Therefore, the nine AF positions are preserved for only 
4% of the transits. 


3.4. Processing of astrometry 


One of the most critical modules of the processing pipeline is 
devoted to process the astrometry, to determine the astrometric 
uncertainties, and to filter the positions that appear to be outliers. 
We describe the different procedures adopted for these tasks be- 
low. 


3.4.1. Coordinate transformations 


We consider here the coordinates associated with a single CCD 
position as produced by IDU for an observed source in the AF 
instrument. These coordinates are expressed in the window ref- 
erence aystem (WRS), and provide the pixel coordinate of the 
centroid of the SSO inside the transmitted pixel window and the 
OBMT reference time of the transit. 


As a first step, the epoch of the crossing of a conventional 
fiducial line on the CCD is computed. 
This corresponds to the exact timing of mid-exposure, and is de- 
pendent on the location of the photocentre of the SSO inside 
the window, on the size of the window itself, the location of the 
window in the focal plane at the time of read-out of its reference 
pixel, and on some more technical aspects such as the binning 
and gating strategy. The OBMT timing of the crossing of the 
fiducial line is effectively the AL coordinate in the WRS. 


The WRS coordinates are then transformed to angular coor- 
dinates in the scanning reference aystem (SRS), whose axes are 
aligned to the AL and AC directions, with the origin at the cen- 
tre of the focal plane (for an overview of these reference systems, 
see Fig. 15 in|Fabricius et al./2016). In this step, the geometric 
calibration of the focal plane 1s applied. In the processing cycle 
producing Gaia DR3, the geometric calibration is among others 
dependent on the source colour, expressed by an effective wave 
number ve. After having performed some tests, it was decided 
that assuming a solar spectrum for all SSOs (corresponding to 
Vert = 0.001561 nm!) was an acceptable approximation. 


A further transformation converts positions from the SRS 
reference system to the centre of mass reference system 
(CoMRS), non-rotating, with the origin at the centre of mass of 
Gaia. A last transformation produces the positions in the BCRS. 
In this step, the relativistic stellar aberration is removed; in other 
terms, the effect of the orbital motion of Gaia is suppressed. The 
outcomes are the position in equatorial coordinates, as seen from 
the centre of mass of Gaia and the associated TCB. 


We stress again that due to the scanning motion of Gaia, the 
TCB of an observation is directly linked to the position of the 
source on the focal plane in the AL direction. The uncertainty 
on the provided TCB is ~1 ms, during which an asteroid moves 
by no more than ~200 pas (a very high value that is only reached 
by some near-Earth Objects). This is negligible with respect to 
the error budget illustrated below. 


3.4.2. Astrometric uncertainties 


We consider a simplified error model that separates uncertainty 
sources that are uncorrelated across a transit from one AF CCD 
to the next (random component) to uncertainties that are con- 
sidered not to vary along a transit (systematic component). This 
scheme is represented by the complete covariance matrix W of 
the transit, 


WwW, 0 oa 0 Ww, W, W, 
0 W> 0 W, W, W, 

W- + , (A) 
0 0 Wo WwW, Ws W, 


where W,, is the covariance matrix of the right ascension and 
declination of the AFn position, which we call the random un- 
certainty of the AFn position, W, is a constant covariance matrix 
throughout the transit, which we call the systematic uncertainty 
of the transit, and 0 is a 2 x 2 matrix of zeros. 

The random component incorporates the uncertainty from 
the centroiding, a term that we call excess noise, and a contri- 
bution from the attitude. Readers familiar with the general astro- 
metric processing by Gaia should note that this is conceptually 
similar to the excess noise defined in AGIS 
[202 1), but its formulation is different. These terms are quadrati- 
cally summed to obtain the total random uncertainty. 

The uncertainty on the centroiding is initially provided as un- 
correlated errors in AL and AC, but after transformation to the 
equatorial coordinate system, the corresponding uncertainties on 
right ascension and declination become highly correlated. How- 
ever, taking the correlation into account, the user can recover the 
precise AL component of the uncertainty. Whereas uncertainties 
in right ascension and declination are typically about 500 mas, 
the real uncertainty in AL is often smaller than 1 mas. 

In AL, uncertainties are at mas level and show the extreme 
precision of Gaia. In AC, the situation is more complex because 
we lack knowledge about a precise position. For sources with 
G > 13, all pixels are binned to a single sample, and the posi- 
tion given corresponds to the centre of the transmitted window. 
The transmitted window is determined by the on-board software 
at the beginning of the transit, so that the object is in one of the 
two central pixels, and is propagated to the next AFs such as to 
keep a non-moving object in the centre. Due to its motion, the 
SSO drifts away from the centre of the window. Therefore, for 
AF CCDs that are reached by the signal later on during the tran- 
sit, the SSO can be anywhere in the window rather than in one 
of the central pixels. By assuming the dispersion of a rectangular 
distribution over the complete transmitted window as value for 
the uncertainty, it is clear that for early AFs, this is an overes- 
timation. However, this approximation of the AC error model is 
not expected to impact the exploitation of the astrometry (e.g. for 
orbit computation) as the AC uncertainties remain ~2—3 orders 
of magnitude larger than the AL uncertainties (further details are 
provided in Sect.|3.4.5). 

For objects G < 13, 2D windows were downlinked to Earth, 
and a 2D centroid fitting was possible. In this case, the uncer- 
tainty in the AC direction is comparable to (but still larger than) 
the uncertainty in the AL direction. 

The second random component to the astrometric uncer- 
tainty, the excess noise, was determined by analysing the orbit 
post-fit residuals of asteroids. It essentially affects bright objects 
(G < 12), so that its origin could be linked to partially resolved 
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Fig. 6: Error model in the AL direction for the SSO astrometry 
in Gaia DR3 as a function of the G magnitude. The total error is 
represented as given by the squared sum of the random and the 
systematic component. The colour represents the data density 
(yellow, lighter colour: higher density). The thick line and the 
two thin lines on each side are the quantiles corresponding to the 
mean and the 1-o level. 


shape effects (Sect.[5). For Gaia DR3, this additional uncertainty 
is given by 

Eex,at, = 0.72 x e0-3578G-10) mas, in AL, 5 
eae = 1.5 x e~932673(G-10) mas, in AC, 2) 
where G is the preliminary magnitude available from the pipeline 
input. 

Finally, the uncertainty on the attitude associated with the 
AGIS solution is derived by analysing the dispersion of epoch 
positions in a magnitude range in which centroiding errors are 
negligible. It contributes both to the random and systematic er- 
rors. The values in Table [3] are the default ones, corresponding 
to regular period of operations (occasional time intervals with 
larger uncertainties are present). 

The resulting uncertainties for all the observations published 
in Gaia DR3 are distributed as illustrated in Fig. [6]as a function 
of the G magnitude for the AL direction. The dominating excess 
noise in the branch of G<12 and the average uncertainties of 
always <1 mas for 10<G<18 and <10 mas at G~21 are striking. 
We show below in more detail how the post-fit residuals to orbits 
represent this quality of the data. 


Table 3: Default uncertainties from the AGIS attitude (mas). 


AL AC 
Contribution to the systematic uncertainty 0.051 2.3 
Contribution to the random uncertainty 0.20 = =1.37 


3.4.3. Filtering of positions 


Filtering the outliers, which are most probably observations un- 
related to SSOs, or data affected by anomalies, is an important 
and delicate task that has been applied both at the level of indi- 
vidual positions (individual CCDs) and at the level of complete 


Article number, page 10 of 34 


transits. The filter definition has required several iterations to op- 
timize the rejection parameters and obtain the cleanest data set, 
minimising at the same time the number of rejected genuinely 
good positions. The whole filtering procedure is described in 
detail in the documentation of Gaia DR3, therefore we do not 
repeat it here. We just recall the broad test categories that have 
been implemented. 


— Test for anomalies at CCD level: For example, samples elim- 
inated or set to zero, attitude of poor quality, close proxim- 
ity of stars, or missing data. Two important rejections occur 
at this stage: all positions by the Sky Mapper (SM) CCD 
columns that are of lower quality are rejected, and all posi- 
tions for objects fainter than G = 22 beyond a reasonable 
range for accurate astrometry. 

— Test for the position uncertainty, within predefined limits 
(function of magnitude). This is based on the study of the er- 
ror distribution (following the model illustrated above) and 
on the identification of clear outliers. 

— Test on residuals with respect to the fit of a linear motion 
during a transit. 


The two dominating rejection reasons (8.6% of the positions) 
are the lack of attitude data and unrealistically large or small un- 
certainties. All the criteria above are responsible for the elimina- 
tion of a fraction of positions that is considered as outliers. An 
extensive analysis has shown that all transits with >2 remaining 
positions can be considered with high confidence as associated 
with a real SSO. Transits that at the end of the filtering have a 
single position are rejected. At the end of the filtering process, 
10.9% of the positions are rejected. 

The astrometric processing is agnostic of the association of 
transits in bundles. It operates only at transit level and treats all 
transits individually and independently. 

As in Gaia DR2, the asteroid astrometry at CCD level of 
Gaia DR3 is provided in the table gaiadr3.sso_observation 
of the Gaia archive. It contains all the data required for its ex- 
ploitation. 


3.4.4. Difference in the processing of astrometry with respect 
to DR2 


A single important difference should be considered by the data 
users that intend to exploit the astrometry at its full accuracy. In 
Gaia DR2, the relativistic light bending was eliminated from the 
observations by applying a correction that assumed the source 
to lie at an infinite distance. After consulting potential users, we 
decided not to apply any light bending correction in Gaia DR3 
and let the users include its effect in the models they use in their 
orbital fitting programs. Depending on the solar elongation, the 
difference of positions due to light bending can be about a few 
milliarseconds, which means that it is relevant to exploit the full 
accuracy. 

Two others choices were made that impact the data volume, 
but not the usage of data. First, bright SSOs have residuals higher 
than the expectation based on formal uncertainties (see Eq. (4). 
In Gaia DR2, we took the simple approach to exclude all SSOs 
with G<10. For Gaia DR3, the higher uncertainty is introduced 
in the error model by the excess noise (Sect. |3.4.2), and bright 
objects are preserved in the output. Finally, a small number of 
sample windows are truncated. It is not straightforward to pro- 
vide a quantitative estimate because 0.5% of the transits has a 
truncation flag set. Of these, only ~10% are probably really trun- 
cated, however. Their positions were discarded in Gaia DR2, but 
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Fig. 7: General illustration of the motion of an asteroid on the focal plane (left panel) and how this is translated into a position 
by the on-board windowing strategy. The horizontal axis represents the AC direction, and the vertical coloured bars are the pixels 
in the AC direction. The vertical axis represents the different CCD strips (their separation is conventional in this scheme). The 
coloured squares connected on the black line represent the real positions of the asteroid in the different CCD strips of the SM and 
AF instrument. The rectangular coloured boxes encompassing 12 AC pixels represent the assigned windows. The dashed black line 
represents the motion that a hypothetical star would have, starting from the same position as the asteroid in SM. The assignment of 
the windows in AF is such that their AC positions, corresponding to their centres (coloured circles), closely follow the dashed line. 
This AC position is transmitted to the ground. The horizontal light grey lines represent the uncertainties on the positions in AC, 
which correspond to the dispersion of a rectangular distribution over the assigned window. The right panel shows how this translates 


into positions on the sky plane. 


are preserved in Gaia DR3 because a dedicated study has shown 
that the quality of their astrometry is not degraded. 


3.4.5. Interpretation of positions during a transit 


The users of asteroid astrometry will see that the positions pro- 
vided along a transit for an asteroid in general follow a zigzag 
pattern on the sky whose average displacement is related, but 
does not correspond exactly, to the asteroid proper motion. In 
the scheme of Fig.{7|(left panel), we explain this effect. It results 
from the on-board windowing strategy. 


When a source enters the focal plane and is detected in a Sky 
Mapper CCD (SM), a window is assigned to the source such that 
it is centred in the window within one pixel. In the case of a star, 
its image drifts across the subsequent window strips due to the 
satellite spin (dashed line). Due to precession, it will also drift in 
the AC direction by a few dozen pixels over the entire transit. Its 
trajectory is therefore not exactly vertical in the plot. 


The windows assigned in the SM strips are propagated 
through the AF strips assuming that the motion to be tracked cor- 
responds to that of a star. Since window shifts can only be per- 
formed by integer pixels, the shifts from one strip to the next will 
alternate between two adjacent integer values, causing a kind of 
zigzag motion of the assigned windows in AC direction. 


In general (i.e. for G > 13), windows are fully binned in AC 
direction, and the only astrometric information we have in AC 
direction is the window position and that the object is somewhere 
inside the window. Thus, the best estimate of the position in AC 
is the window centre, causing the derived positions to follow the 
same zigzag motion as the window. The uncertainty on the po- 
sition is consequently given as the dispersion of a rectangular 
distribution over the assigned window (shown as the horizontal 
grey bar in the figure). 


As the AC motion of the real asteroid is totally independent 
of the window propagation, its signal can be severely truncated 
and the asteroid can leave the window. The derived position is 
thus affected by a large bias. In the scheme, this is represented 
by the fading of the circles when the asteroid reaches the window 


edges. AF7-9 positions will probably be rejected as unreliable by 
the astrometric processing. 


After they are transformed onto the sky (right panel), the 
different windows nearly overlap in AL, but show shifts in AC 
aligned on integer pixels. In the figure, a small AL motion (com- 
ing from the asteroid proper motion) is present in addition to 
the zigzag motion. It should be emphasised that after rotation to 
RA-Dec, the zigzag motion may visually appear in both coor- 
dinates, and that the uncertainties will seem to be large both in 
right ascension and declination, and thus (only apparently) mask 
the extreme precision of the position, which is still present in the 
AL direction. 


While it can be important to understand the origin of this pe- 
culiar pattern, its influence on practical applications (e.g. orbit 
fitting) is negligible, as the spurious AC motion and its fluctu- 
ations will be zero within the given uncertainties in AC. Con- 
versely, the AL position is very accurate and also represents the 
proper motion of the asteroid in that direction with high preci- 
sion. 


3.5. Orbit computation 


The computation of orbits from the Gaia astrometry alone has 
not been run with the core of the pipeline, but as a post- 
processing task. In the Gaia DR3 archive, orbits are published 
in a specific table named gaiadr3.sso_orbits. 


Orbital fitting is performed based on Gaia-centric astromet- 
ric data in right ascension and declination at CCD level, using the 
corresponding error model with non-diagonal covariance matrix 
as weights, as derived in Sect. [3.4.2| Only known asteroids are 
considered. 


For this task, we adopt the usual procedure consisting in de- 
termining the corrections to the initial orbital elements by solv- 
ing a linear system of equations for each asteroid, while min- 
imising the residuals between the computed and the observed 
positions. To determine the computed values as well as the re- 
quired partial derivatives, the equation of motion for each source, 
including the relativistic contribution, is integrated simultane- 
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Fig. 8: Computed uncertainty of the semi-major axis o-, for the orbits published with Gaia DR3 as a function of the semi-major axis 
itself. Colour represents the density of objects computed in the bins (600x500 bins in the axis range). 


ously with the variational equations 
[1969p 

The positions and velocities of the planets of the Solar Sys- 
tem are obtained from the highly accurate dynamical model IN- 
POP10e (consistently with the whole Gaia software framework), 
which includes the eight planets, the dwarf planet Pluto, and a se- 
lection of 343 asteroids Deram et al.[2022). 
The heliocentric positions of the asteroids are computed from 
numerical integration including perturbations from eight planets 
and Pluto; relativistic corrections are made with the parametrised 
post-Newtonian (PPN) formulation. 

Outlier rejection at the observation level has been imple- 
mented as in (2003). We included a test on tol- 
erance and on the maximum number of iterations. 

The reference epoch for the orbit solution can in principle 
be chosen arbitrarily and could be the same for all asteroids; 
however, a more precise result is obtained with a reference epoch 
halfway on the observational arc. This value is provided in the 
Gaia archive. 

The final output of the procedure is the new improved state 
vector in the ICRF3 system at the reference epoch, together with 
its corresponding covariance uncertainty matrix. These quanti- 
ties are also transformed into heliocentric elliptical orbital el- 
ements with the associated covariance matrix at the reference 
epoch (see Sect. [2.2] for the details about the reference systems 
used). 

Residuals are expressed in right ascension and declination, 
but with correlation due to the orientation of the Gaia scan- 
ning plane during the measure. When the scanning direction on 
the sky is known, these residuals can also be expressed in the 
independent along-scan and across-scan directions, that is, in 
the (AL, AC) plane, without a correlation. General validation 
is performed on the residuals expressed in the (AL, AC) coor- 
dinates, and more particularly, in the most precise AL direction 
(see Sect.|5.1p. 

The orbital fit implementation was validated by various com- 
parisons of the state vector, orbital elements, and their associ- 
ated covariance matrices. This includes internal and external data 
and codes, in particular, a comparison with orbits derived inde- 
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pendently from the preliminary astrometry provided by IDT, the 
astorb and JPL databases, covariance, and the NIMA code (J. 
Desmars, priv. com.) for orbit determination and propagation. A 
final validation into the DPAC frame has been also performed by 
using the orbfit code and is discussed in[Babusiaux| (2022). 

The number of known asteroids for which the orbit fit was 
run is 156 801. The procedure did not converge for 451 sources. 
The final number of orbits is 156 350. 

Finally, not all orbit solutions are meaningful because only 
a small number of astrometric data is available and the cov- 
ered time span is limited. We decided to discard extremely un- 
certain orbits, defined as being based on <20 observations, or 
those for which the astrometry covered <60 days. Another 29 
orbits were discarded by a direct threshold on relative uncer- 
tainty 7,/a > 9 10~*. A total of 1609 orbits were finally rejected 
(~1%). The final number of orbits is 154741. Figure[8]shows the 
general distribution of the uncertainty on the semi-major axis for 
the published orbits. 

A total of 23 327 388 observations were processed, grouped 
as 3212676 transits. Out of this grand total of observations, 
23 031703 were not rejected by the fitting process. The rejected 
observations (~1.27% ) mostly belong to the 1609 eliminated 
orbits. The set of retained orbits contains 23 167 198 observa- 
tions grouped into 3 192 098 transits, so that the proportion of 
rejected observations from this set of objects corresponds to only 
~0.58%. 


3.6. Processing of photometry 


G-band photometry is obtained from the signal measured by the 
unfiltered AF instrument, the as was same used for astrometry. 
For the SSOs, the G-band photometry is not derived in the frame 
of their specific processing modules, but similar as in Gaia DR2, 
by an upstream system (PhotPipe) that treats and calibrates pho- 
tometry for all the sources observed by Gaia, producing fluxes 
and their errors for all the AF CCDs (Riello et al.J2021). 

In the frame of the SSO processing, only the weighted aver- 
age of the fluxes at transit level is computed. In this process, 
measurements that are rejected by the astrometric processing 
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(mentioned above) are also rejected. The rationale for this is the 
fact that problems with the signal that prevent the computation of 
an accurate position (most often truncation, attitude problems, or 
other anomalies) also produce problematic flux measurements. 
Other specific flags set by upstream systems are also reasons for 
flux rejections. The details of these rejection criteria are given in 
the Gaia DR3 

The average G-band flux per transit was computed using the 
weighted average, where weights correspond to the inverse of 
each flux variance. The quadratic sum of the flux errors per CCD 
was used to obtain the flux error for the transit. The average G 
-magnitude value was computed by using the average G-band 
flux and the G -magnitude zero-point provided by the general 
photometric calibration. 

The flux, its uncertainty, and the magnitude can be accessed 
in the table gaiadr3.sso_observation in the fields g_flux, 
g_flux_error and g_mag. While it is recommended to use the 
flux uncertainty, conversion into error on magnitude can be per- 
formed by taking the zero-point error provided in Sect. 5.4.1.3.2 
of the DR3 online documentation into account. The average 
value repeats identically for all observations in the transit. As 
a post-processing task, all magnitudes G>21 or with errors >0.2 
were considered as unreliable and were eliminated. 


4. Evaluation of the sample of unmatched sources 


As explained in Sect.|3.1| the search for unmatched sources was 
performed with orbital elements as known at the end of 2017. 
Meanwhile, many new orbital elements have been computed and 
published. In early February 2022, we did a search to identify 
the unmatched sources using the most recent orbital elements 
published by the Minor Planet Center. 

Because the list of unmatched sources contains 3541 tran- 
sits, and the file with orbits contains about 1.2 x 10° orbits, this 
results in a total of ~4 x 10° residuals that were to be confirmed, 
each requiring a numerical integration of the orbit with full per- 
turbations by all major planets. To reduce the number of com- 
putations to a manageable quantity, some optimisations were es- 
sential. 

Therefore, we divided the sky into tiles of roughly one 
square degree each. In a first iteration, we computed geocentric 
ephemerides of all asteroids with a published orbit at one-day in- 
tervals, without perturbations, to determine the tile in which the 
asteroid is predicted to be for each Gaia transit of an unmatched 
source. Correspondence to the transit positions was tested in a 
tile set that included the central tile plus surrounding tiles (with 
at least one corner in common). This compensates for the paral- 
lax effect between Gaia and the geocentre, and for the absence 
of perturbations. Potential identifications are all asteroids found 
in the tile set at the epoch closest to the epoch of the transit. Far 
from the ecliptic, this can result in fewer than ten candidate iden- 
tifications, which grows to >1000 identifications in the ecliptic 
region. 

The second step was to compute the residuals between the 
Gaia position and the position from ephemerides for each can- 
didate identification, but now Gaia-centric, computed with full 
perturbations, and for the exact epoch. Since the motion of an 
asteroid in the course of a transit is small compared to the uncer- 
tainty in position in AC, we used only one position per tran- 
sit. We applied this procedure three times, with three differ- 
ent values of the threshold for a tile to be a neighbour: 2 de- 
grees, 5 degrees, and 10 degrees. In each iteration, only the still- 
unidentified sources were considered. The last iteration did not 
give any additional identification. 


Next, we defined some criteria to accept an identification. 
The key point in this is to rotate the residual in right ascension 
and declination to the directions parallel and perpendicular to the 
line of variation. The line of variation is defined as the line on the 
sky that results from varying the mean anomaly of the asteroid 
without changing the epoch or the other orbital elements. The ra- 
tionale behind this is that errors in the orbital elements will result 
in periodic errors in the position of the object, except for an error 
in the semi-major axis (or the period), which will result in a drift 
in mean anomaly, linear with time. This will cause the error in 
mean anomaly at the epoch to have the largest impact in the error 
in position. To find suitable thresholds, we performed a similar 
exercise using an orbital elements file from 2014, in which some 
of the now numbered minor planets still had rather poor orbits. 
Thus, we were able to determine how residuals behave for identi- 
fications that are known to be correct compared to incorrect iden- 
tifications. For incorrect identifications, we found differences in 
motion that were evenly distributed between 0 and more than 
500” per day, while for correct identifications, we found that the 
differences rarely exceeded 2” per day. For correct identifica- 
tions, however, positions may be off by up to 10’ along the line 
of variation. Perpendicular to the line of variation, residuals do 
not exceed 10”. No incorrect identification had residuals smaller 
than 10” perpendicular to the line of variation, with a difference 
in speed smaller than 2” per day at the same time. Thus, it is pos- 
sible to distinguish between correct and incorrect identifications 
with only a very small risk of an incorrect identifications if the 
speed can be measured, that is, if there are at least two transits. 

For unmatched sources with at least two transits, we there- 
fore set up as limits 400” for the residual along the line of vari- 
ation, 20” for the residual perpendicular to the line of variation, 
and 7”/day for both components of the difference in motion. This 
last rather high value is justified by the fact that if an object has 
only two transits that occur less than 0.25 day apart, for instance, 
the uncertainty on the motion is larger than the expected differ- 
ences in speed. For unmatched sources with only one transit, we 
had to be more severe on the residuals on the positions because 
the speed could not be computed. We set 5” as limit along the 
line of variation and 1.5” perpendicular to the line of variation. 

With these criteria, we found an identification for 712 out 
of the 1320 unmatched sources. Some of the unmatched sources 
turned out to be the same asteroid, so that in total, they repre- 
sent only 567 different asteroids. Strangely, the highest residual 
we found was only 13”, rather than the expected several hundred 
arcseconds. This shows the lack of preliminary orbits in the 2022 
Minor Planet Center orbital elements files. The derived identifi- 
cations are made available in the auxiliary data web page of Gaia 
DR3. 


5. Astrometric performance 
5.1. Orbit quality 


The large differences in both data volume and time coverage for 
each asteroid source means that the quality of orbits computed 
from Gaia DR3 astrometry varies from excellent to very approx- 
imate. This is shown by the values of the semi-major axis uncer- 
tainty o, in Fig. [8] 

A bimodal distribution of the uncertainty appears. This fea- 
ture is more clearly visible from the frequency distribution of 
the relative uncertainty (Fig.|9), in a limited range of semi-major 
axis (here from 2.2 to 2.5 au). This signature can also be detected 
in similar distributions obtained with ground-based observations 
(Desmars et al.|2013), reflecting a combination of short observa- 
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Fig. 9: Distribution of the relative uncertainty on the semi-major 
axis for the computed orbits in the range of a between 2.2 and 
2.5 au. The right peak of high values in general contains orbits 
obtained from a small number of observations over a short ob- 
servational arc. 


tion arcs and small astrometric data sets. This can be confirmed 
by the statistics derived for the two populations, assuming that 
they are separated at the threshold value 7,/a = 10°°. The 
average number of non-rejected observations is 147.5 and 79.1 
for the lowest and highest uncertainty, respectively, while the arc 
lengths are 771 and 380 days. 

An interesting question related to the quality of the orbital 
solution from astrometry in Gaia DR3 is how it compares to the 
orbits obtained from larger data volume and longer time-span, 
including the whole set of astrometric measurements available 
from the ground. In particular, we recall here that the time span 
of Gaia DR3 is just 34 months, which corresponds to the or- 
bital period of an asteroid with a semi-major axis of slightly less 
than 2 au. To address this question, we have obtained the most 
recent orbital solution (updated on 7 March 2022) from the JPL 
database of orbits“]for all objects for which we have an orbit. We 
also used the related absolute magnitude H provided by the same 
source to categorise the comparison with respect to the absolute 
brightness. 

In Fig.[10|we plot the ratio of 7, from JPL and the Gaia DR3 
value as a function of semi-major axis and H. A general trend of 
increasing accuracy with brightness appears for MBAs. While 
most of the orbits lie below the line of equivalent accuracy, as 
expected from the long time coverage and data value for the gen- 
eral population, it is interesting to note that 8736 asteroids reach 
a better accuracy when Gaia-only data are used, including some 
NEOs, several Jupiter Trojans, and TNOs. 

The post-fit residuals to the orbit adjustment are a final sen- 
sitive test for the quality of the astrometry. Figure [11] presents 
the global view of the residual dispersion. The distribution is 
strongly dispersed in the AC direction, as expected, due to the 
low accuracy of astrometry in the across-scan direction with re- 
spect to the AL direction. 

Extreme values of residuals correspond to the less accu- 
rate orbit of faint objects with fewer observations. Although the 
range of values appears large, the core of the distribution is very 


4 https://ssd. jpl.nasa.gov/tools/sbdb_query 
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Fig. 10: Ratio of 0, from JPL and from Gaia DR3 as a func- 
tion of semi-major axis (top panel) and of absolute magnitude 
H (bottom). Colour represents H (top) and semi-major axis (bot- 
tom). The horizontal line shows equal uncertainties (ratio=1). 


compact. This is visible in the histograms of Figure [11] where 
the central peak appears. The global histogram, however, does 
not convey complete information on the astrometric accuracy. To 
highlight the quality of Gaia astrometry as a function of bright- 


ness, we chose the same approach as was used in|Gaia Collabo- 


(2018), and we represent the dispersion (represented 
by the standard deviation) of the residuals over a transit (Fig. [12). 


The average and quantile lines clearly illustrate the trend. 
The average value reaches a bottom plateau at the exceptional 
value of ~0.25 mas over a range of three magnitudes, from 
G = 12 to 15 mag. Up to G = 17 mag, the dispersion re- 
mains at submilliarcsecond level. This is in general agreement 
with Fig. (6| although the orbital post-fit residuals can still re- 
flect some systematic effects that are not yet taken into account 
in the modelling (see Sect. [5.2.2). In this range, the accuracy of 
Gaia DR3 is clearly higher than that of Gaia DR2. Moreover, 
the transitions in the error values that were visible in Gaia DR2 
at G = 13 and G = 16 mag have now disappeared as a result of 
the increased global quality of the calibration. Figure[13]|shows 
that the improvement reaches a factor of about two and appears 
fully consistent with the robust estimate of the standard devia- 
tion for the general astrometric processing of stars in Fig. A.1 of 


Lindegren et al.|(2021). 
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Fig. 11: Distribution of the post-fit astrometric residuals in the 
AL,AC plane. The vertical lines mark the 1-sigma quantiles 
computed for bins of 90 mas on the AC axis. The frequency dis- 
tributions along the two axes are reproduced in the middle and 
bottom panel (in AL and AC, respectively). The corresponding 
standard deviations are 5.15 mas and 270.14 mas. 


5.2. Shape and size effects 
5.2.1. Binary systems 


As foreseen (Tanga et al.!2008}|Pravec & Scheirich}2012}/Tanga 
& Hestroffer)/2012), one of the most interesting applications ex- 
pected for the accurate astrometry by Gaia is the possibility of 
detecting satellites of asteroids. As the orbital fit tends to con- 
verge to the trajectory followed by the centre of mass of the ob- 
ject, residuals can contain the signature of an asymmetric light 
distribution between the primary and the secondary component 
of an unresolved binary. The expected amplitude and periodicity 
of these residuals can fall in a range that is accessible to Gaia, 
and they can also cover a range of component ratios and sizes 
that is not accessible to other commonly used techniques (in- 
cluding radar ranging, adaptive-optics imaging, and photometry 
of mutual eclipses or occultations; 2015). 
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Fig. 12: Dispersion of post-fit residuals in AL for each transit. 
The colour represents the data density. The deep blue line shows 
the mean, and the light blue lines show the quantiles correspond- 
ing to 1-sigma, computed over 50 bins in the interval of G mag- 
nitude from 8 to 21 mag. 
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Fig. 13: Average and 1-o quantiles for the AL dispersion of the 
post-fit residuals for all transits in common between DR2 and 
DR3. In contrast to Fig. the single data points are not repre- 
sented. The background curve (light blue) represents the distri- 
bution computed for Gaia DR2. 


The amplitude of the wobbling w, that is, the maximum dis- 
tance between the system photocentre and barycentre seen by 
an observer, can be easily estimated by assuming spherical com- 
ponents of identical albedo and bulk density, characterised by a 
diameter ratio k = D>/D,, hence a mass ratio g = k° 
ifer et al.[2010). With this notation, the ratio of their illuminated 


surfaces is proportional to g?/* and determine the position of the 
photocentre. It is then easy to show that 


1 1 
l+g?? 1+q'! 


a, (6) 


where a is the semi-major axis of the mutual orbit. In this ap- 
proximation, no phase effects are introduced so that the individ- 
ual photocentre offsets for each component are neglected. 

A tempting opportunity to look for binary-related photocen- 
tre wobbling in Gaia data came from the discovery of the bina- 
rity of the gus (4337) Arecibo through two lucky occultation 
events (Fig.[14) 20.78 days apart in May and June 2021 (Gault| 
fet al.||2022). The two occultation chords for each event pro- 
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Fig. 14: Size and relative position of the primary and secondary 
components of (4337) Arecibo as derived from the occultation 
(in blue) of 19 May 2021 Cleft), and 9 June 2021 (right) on the 
plane of the sky in the equatorial reference. The red crosses and 
the orbit are derived from the model described in the text. 


vided diameter estimations from a circular fit of the components 
(D,=24.4+0.6 km, D2=13.0+1.5 km; k,.¢=0.55) and two accu- 
rate astrometric measurements of their relative positions, with an 
apparent separation of the components of 25 and 32 mas. In the 
absence of any additional information, it is not possible to de- 
rive information about the mutual orbit from these two relative 
positions. At the given apparent separation, however, it is clear 
that the object does not appear as resolved to Gaia. The signal 
of the two components falls within one pixel (60 mas AL). How- 
ever, given the extreme centroiding accuracy in the AL direc- 
tions, perturbations corresponding to the photocentre wobbling 
can be expected. 


Gaia DR3 contains astrometric data for (4337) Arecibo in 
38 transits. The G magnitude is present for 36 transits, so that 
a photometric inversion was attempted with the genetic algo- 
rithm used in (2019). This provided a period 
of 32.972823 hours and pole ecliptic coordinates (A, 8)=(271°, 
68°). Given the long time span covered by the observations (881 
days), the period appears rather well constrained (at a level prob- 
ably better than 1s), while the typically estimated uncertainty 
of the pole direction is ~10°. Moreover, the rotation period is 
compatible with a preliminary unpublished light curve obtained 
from ground-based observers in the weeks following the occul- 
tation (32.85+0.38 hours; Behrend, R. et al., private communica- 
tion). We stress here that we do not expect to directly find indica- 
tions of the binary nature of the asteroid in the sparse photometry 
alone. 


In the Gaia data set, we looked for the longest possible se- 
quence of consecutive transits and found one composed of 13 
transits over 2.3 days, 12 of which are consecutive (correspond- 
ing to six rotations of Gaia). During this short time span, the 
orientation of the satellite orbit with respect to the observer does 
not change appreciably. This sequence is therefore a good can- 
didate to search for any periodicity in the residuals that could 
suggest the presence of wobbling. In addition, during a single 
transit, the astrometry is expected to be affected by a systematic 
displacement, if present, by the same amount. It is thus possible 
to exploit the average of all residuals collected over a single tran- 
sit for a more robust estimate. We also assume that the standard 
deviation of the residuals represents the typical uncertainty. 

As the AL direction conveys the accurate astrometry, only 
AL residuals were taken into account. Therefore, any mea- 
sured wobbling would be caused by the component of the two- 
dimensional photocentre-barycentre shift (as it appears to the ob- 
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server, on the sky plane) in the direction of the scan (AL) for each 
transit. 

The analysis of the average residuals as a function of time 
within the selected sequence clearly suggests a systematic fluc- 
tuation. Figure |15|shows that fluctuations occur on a timescale 
compatible with the rotation period. Based on the considerations 
above, the compactness of the system is suggestive of a satellite 
revolution period that is synchronous with the primary rotation, 
which would explain the compatibility between the astrometric 
and the photometric periods. 

With the elements above, we can now show that Gaia is able 
to provide an orbital solution for the system. We assumed that 
the most robust parameter a priori available for the model is the 
rotation period, and based our analysis on the relative astrometry 
derived from the occultation and the AL transit residuals, with 
their error bars. Only the occultation chords, but not the derived 
relative astrometry that we exploit, are sensitive to the absolute 
sizes D,.. For the same reason, while the wobbling recorded by 
Gaia in the residuals depends on the size ratio k through Eq. [6] 
the occultation results do not. 

To model the wobbling, we exploited Eq. [6]and the position 
angle of the scan, available with the astrometric data, to com- 
pute its component projected in AL. The model also requires the 
ephemeris of the asteroid as ancillary data to correctly repro- 
duce the observation geometry, including the light—time delay to 
be considered for each observation epoch. 

We defined a target function as the squared sum of the O- 
C of the model-derived astrometry with respect to the measured 
astrometry. We considered the pole coordinates (A, 6), the semi- 
major axis a, and the initial rotation phase at an arbitrary refer- 
ence epoch (set to the first occultation event) as free parameters 
of the model. A downhill simplex algorithm (Nelder-Mead) was 
run for its minimisation. Our result provides the values 


a= 49.9 + 1.0 km, (A,B) = (261° + 3°, 60° + 3°). 


The pole coordinates remain compatible with those derived 
by photometric inversion that were used as initial conditions. 
The fit to the occultation data is strikingly good (Fig. (14), al- 
though the final result shows a small discrepancy in the occul- 
tation astrometry, especially for the first event. However, the 
chords of the first event are very similar to each other, and er- 
rors related to shape assumptions (spherical shape hypothesis) 
could show up at the milliarcsecond level. 

The wobbling amplitude following Eq. [6] could reach 8.5% 
of the object separation, or about 2.7 mas at the distance of Gaia 
(2.24 au) for the observations exploited here. This amount is 
44% of the apparent radius of the primary component. It is there- 
fore much larger than a photocentre shift due to shape effects 
alone for the phase angle at the same epoch (14.3°). 

The measured wobbling is strongly reduced by the projection 
in the AL direction (Fig. [16). However, our fit attempts show that 
it cannot be fitted unless a value of k is assumed that is lower by 
~ 35% with respect to the nominal value derived by occulta- 
tions (ky, = 0.35 ko.-). With this correction, the model appears to 
reproduce the observations pretty well (Fig. [15] bottom panel). 

The apparent contradiction between k,-- and ky, can have 
many origins that are all equally interesting, implying that the 
mass of the companion is lower than what is estimated from the 
geometric size ratio provided by the occultations. A first possi- 
bility is a non-spherical shape of the components with a flatten- 
ing that is more pronounced in the case of the satellite. In this 
case, because the occultation chords for the second event appear 
to constrain mainly the equatorial radius of the two components, 
koce WOuld not represent the volume ratio well (the first event 
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Fig. 15: Residuals to the orbital fit of (4337) Arecibo (blue dots). They are obtained from the average of single-observation residuals 
over each transit. The error bars are given by their standard deviation. In the top panel, the dashed grey line is not a model fit, but 
a simple overplotting of a sinusoid of the period derived by photometry, adjusted in amplitude (0.8 mas) and phase to the data. In 
the bottom panel, the same data are shown with the residuals predicted by the optimised binary model described in the text (red 


crosses). 
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Fig. 16: Orientation of the modelled orbit, projected on the equa- 
torial reference (RA* indicates that the factor cos(Dec) is in- 
cluded), at the mid-epoch of the sequence of observations ex- 
ploited to model the binary (4337) Arecibo. The positions of the 
satellite at the epoch of each observation (red dots) are shown, 
together with their projection on the instantaneous direction of 
the scan (AL, red lines). The measured photocentre wobbling is 
proportional to the component of the binary separation in the AL 
direction. 


does not constrain the size ratio well in any case). A second pos- 
sibility, relevant for the formation mechanisms, is that the com- 
panion could have a lower bulk density. In either case, the mass 
ratio g would be reduced. 

Of the two options, shape flattening is likely to be favoured 
because two additional other elements support it. First, the ther- 


mal infrared diameter (Mainzer et al./2011) is 32% smaller than 


the surface-equivalent diameter from the occultations. There- 
fore, the equatorial size constrained by the occultations would 
rather be a maximum shape extent, while the polar radius could 
be smaller. Second, when we apply Kepler’s law to our best- 
fitting parameters, the density is too low when the nominal oc- 
cultation diameter is used. The size must be reduced by 19% to 
increase the density to a minimum value 1 g/cm?. A low density 
like this is expected for this object, which belongs to the Themis 
family, whose potentially water-rich, highly porous members are 
expected to have a density below 1.3 g/cm? (Marsset et al [2016). 

Future observations, by Gaia and from the ground, should 
be able to better characterise this system and also investigate the 


possible role of surface scattering and photocentre shift of each 
of the components. 


5.2.2. Photocentre wobble for (21) Lutetia 


The decrease in accuracy of Gaia astrometry towards the bright 
magnitudes (also seen in Gaia DR2 
2018) is suggestive of the fact that shape and size effects could 
affect the centroid determination and degrade the performance. 
This effect was also present in the astrometry of the Gaia pre- 
cursor HIPPARCOS (Hestroffer} 1998). 

While a detailed and systematic search of this effect in Gaia 
DR3 asteroid data is beyond the scope of this article, a single 
case can already provide interesting indications and can in par- 
ticular show whether a discrepancy between the photocentre and 
the centre of mass of the asteroid is detectable in the astromet- 
ric residuals with respect to the fitted orbit. Asteroid (21) Lute- 
tia with its detailed shape model (sect. 6.2.1) is especially well 
suited for this study. Lutetia can be approximated as an ellip- 
soidal asteroid with total axial dimensions of (12110175) km 
and a mean radial distance from the geometric centre of 49 km. 
By using the high-resolution non-convex shape model resam- 
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Fig. 17: Photocentre and barycentre positions of asteroid (21) Lutetia (brown and blue bullets) as projected in the direction of Gaia 
(normal to the image plane) for three epochs. The z-axis of the equatorial coordinate system is denoted by the green bars. North 
points to the top. The Sun illuminates the asteroid from the directions denoted by the red bars. From left to right, the solar phase 
angle takes the values of 16.7, 24.2, and 20.9 degrees, and the photocentre-barycentre offset takes the values of 7.42, 5.96, and 


10.91 km. 
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Fig. 18: Comparison between the average transit residuals for (21) Lutetia in the AL direction (x-axis) and the computed photocentre 
shift from the light scattering model (y-axis; left). The error bars represent the standard deviation of the residual across the transit 
(often smaller than the symbol size). The right panel shows the equivalent plot after correcting for the for the computed shift. The 
lines represent the displacement of each data point with respect to the left panel. 


pled at five-degree resolution, we computed the photocentre- 
barycentre offset for the 23 epochs of the GDR3 photometry (for 


the computation of the photocentre, see for example 
2015). 

Example offsets are illustrated in Fig.{17|for the illumination 
and observation geometries in three representative epochs, corre- 
sponding to phase angles of 6.7, 24.2, and 20.9 degrees. The cor- 
responding photocentre offsets were 7.42, 5.96, and 10.91 km. 
At the distance of Gaia at the epoch of the observations (3.15, 
2.45, and 2.78 au, respectively), the angular offsets are 3.4, 3.3, 
and 5.4 mas, respectively. 

We exploited the position angle of the scan to project the 
computed offset of each epoch in the AL direction, deriving a 
prediction for the offset that Gaia should have measured. We 
then compared it to the orbital post-fit residuals of Gaia obser- 
vations, always in the AL direction. In this case, the residuals 
are represented by the average of the residuals of the individual 
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observations of each transit. The orbital fit was obtained by the 
same approach as used in Sect. [5.3]and included all the astrome- 
try available at MPC. 


The result is illustrated in Fig.{18} It shows a clear correlation 
between predictions and observations (left panel), with a certain 
scatter. 


We then subtracted from the Gaia astrometry the computed 
photocentre offset, and then applied the same orbital fitting pro- 
cedure, with the subsequent analysis of AL residuals. As ex- 
pected, now the distribution of residuals is much more compact 
(right panel), and the correlation with the computed shift has 
disappeared. The final residuals are distributed around zero with 
a standard deviation of 1.2 mas. This remaining scatter is rel- 
atively large with respect to the brightness of (21) Lutetia, ob- 
served by Gaia at G~13. It suggests that a margin of improve- 
ment to the photocentre model probably exists. Concavities at 
large phase angles (particularly relevant in the third panel of 
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Fig. (17), where shadowing effects can enhance the offset by a 
large amount, play a relevant role. The correct modelling of the 
light-scattering properties of the surface is probably critical in 
these situations. Albedo variations are also not considered, but 
might contribute to the residuals. 


Despite these limitations, this example provides the com- 
pelling evidence that an asteroid in the 100 km class can exhibit 
a rather large photocentre offset that is strongly dependent on 
shape details and reaches an order of 20% of the average radius 
at the phase angles of the observations by Gaia. At the level of 
sensitivity of Gaia astrometry, it should be possible to recover 
this effect in Gaia DR3 on asteroids that are probably two to 
three times smaller. 


5.2.3. Pluto and Charon 


We present here the peculiar case of an emblematic resolved bi- 
nary system: the Pluto and its major satellite Charon. The Pluto 
system is regularly seen by the Gaia optical system, in the same 
way as any other SSO. Pluto is a bright source of magnitude 
G ~ 14.5 and is easily detected. The observations were pro- 
cessed in the same way as for any other SSO. However, the 
largest satellite of Pluto, Charon, is at a brightness G ~ 16.5 mag 
and separated on the sky by at most 1”’ from Pluto. Gaia can de- 
tect both provided the projected separation on the scan direction 
(AL direction) is larger than at least 0’’25, so that the on-board 
detection could resolve the system into two independent sources. 
Over the time interval of DR3, 17 resolved and 6 unresolved pas- 
sages were expected. As explained inB.2| there is no particular 
procedure in the identification step to flag each of the compo- 
nents of a resolved system. They are matched to the system, and 
all observations go through the pipeline. Only at the very end of 
the astrometric solution can we see how many have survived and 
assign the solutions to either Pluto or Charon. In this particular 
case, the G magnitude allows us to obtain a flawless identifica- 
tion, but the comparison to the computed positions would work 
as well. 


We found that six transits were successfully resolved for 
Pluto and Charon, thus providing 44 accurate absolute positions 
for each body. We plot in Fig. [19]the relative positions of Charon 
as referred to Pluto for these six transits, together with an outline 
of the apparent orbit at the times of the first (a) and last (d) tran- 
sit. Pluto is at the centre of the plot. The (a) and (d) observations 
combine transits from both fields of view, while a single transit 
is available in the (b) and (c) passages. The green triangles are 
only for AFI and were singled out because of a frequent offset 
in the relative positions for these observations (in the AC direc- 
tion, therefore irrelevant for the astrometric accuracy of these 
measurements). The red triangles in (a) are stacked for 13 obser- 
vations and in (d) for 14 observations. There are 6 observations 
for (b) and (c) in the upper triangle, including one AF1. The 
lower (c) triangle is AF8, and it is offset compared to the other 
six in this transit. The arrows indicate the direction of scan (AL) 
and its perpendicular in the AC direction. The successful reso- 
lution of the binary system occurs when Pluto and Charon are 
favourably oriented with respect to AL, that is, when their sepa- 
ration projected on AL is large. We have an excellent accuracy in 
this direction, and as expected, we find almost all the points on 
or very close to the computed orbital path at the relevant epoch, 
knowing that (c) is in September 2015, and (d) is much later, in 
March 2017. 
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Fig. 19: Relative positions of Charon with respect to Pluto from 
Gaia DR3, with the computed apparent orbits at the epoch of 
the first (a) and last (d) observation. The red triangles show the 
positions in AF2-AF9 and the green triangle shows AF1. The 
arrows indicate the AL-AC directions (see text for the details). 


5.3. Yarkovsky effect detection 


The Yarkovsky effect is the most important non-gravitational 


perturbation acting on small Solar System bodies (Vokrouhlicky 


2000). This secular perturbation produces a drift in the 
semi-major axis of the objects, changing the orbit of small aster- 


oids over millions of years. It is now considered to be responsible 
for their migration from the main belt to the near-Earth region, 
and it represents the key to understanding the evolution of as- 


teroid families (Bottke et al.]2001 2015 
2022). 


The Yarkovsky effect is proportional to the inverse of the di- 
ameter (larger on smaller bodies) and it depends on several phys- 
ical quantities, such as the thermal inertia, the Bond albedo, the 
density of the object, and the rotation period, which are usually 
unknown. As a consequence, different methods have been devel- 
oped to directly measure the Yarkovsky effect from the astrom- 


etry (Farnocchia et al.//2013}|Del Vigna et al./2018 
2020). 


These methods can easily lead to false detections, especially 
when the astrometry contains errors that are usually hidden in 
the low quality of the data. To avoid this possibility, a detection 
is usually considered valid if the following conditions are met: 


— The signal-to-noise ratio (S/N) of the Yarkosky measurement 
is higher than 3. 

— The ratio of the expected value and the actual value is lower 
than 2, where the expected value is an approximation of the 
value that we would expect to find for an asteroid of the same 
size as the one we consider. 


For a complete description of the validation methods, we refer 
to Del Vigna et al.]2018). The Yarkovsky 
effect has so far been measured for 234 asteroids (source JPL 
Small-Body Databasd’), and all of them belong to the NEO pop- 
ulation. The explanation is easy: we need accurate orbits and 
small objects, and these two characteristics are usually easier to 
be found in the NEO population. NEO orbits are better studied 
because of the impact hazard, and radar observations can also 


> https://ssd. jpl.nasa.gov/tools/sbdb_lookup.html#/ 
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be performed at their close approaches to the Earth. A combina- 
tion of all these points makes their orbits more accurate and less 
prone to errors. 

As already mentioned at the beginning of the section, the 
Yarkovsky effect is also the main key for understanding the evo- 
lution of asteroid families. Families are generated by past colli- 
sions between asteroids. The orbits of the smaller members of 
the families moved from their initial configuration because they 
were perturbed by the Yarkovsky effect. A measurement of this 
effect gives us the age of the family, which corresponds to the 
time of the initial collision (Spoto et al./2015). The Yarkovsky 
effect has so far never been measured in the main asteroid belt, 
mostly because we lack accurate observations for MBAs. 

The orbits and their uncertainties were obtained as a result of 
the validation procedure. An orbit determination fit, independent 
from the one presented in Sect. was performed with a mod- 
ified version of the OrbFit software? More information about 
the independent validation of Gaia DR3 Solar System observa- 


tions can be found in/Babusiaux| (2022). 
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Fig. 20: Semi-major axis (au) vs log 10 of the semi-major axis 
uncertainty (au) for the objects for which the orbit determination 
procedure converged during the validation process. The colour 
bar represents the diameter estimate for each object from JPL 
SBDB. 


The results in Fig. [20] were obtained using Gaia DR3 ob- 
servations alone. They show the semi-major axis and its cor- 
responding uncertainty o,. The colour bar represents the esti- 
mate of the diameter for each object from the JPL Small-Body 
Database (SBDB). o, represents a measure of the quality of the 
orbit. It is clear that the orbits of larger objects are better con- 
strained in the main belt, while it is easier to find very good 
orbits among smaller NEOs, even using Gaia DR3 alone. The 
34 months of observations covered by Gaia DR3 are still a too 
short time interval to detect the Yarkovsky effect for MBAs. We 
do not expect to be able to use Gaia DR3 observations alone, but 
itis worth noting that some objects still reach extremely small 
orbital uncertainties, as shown in Fig. 

To fully exploit the data, we combined Gaia DR3 ob- 
servations with all the available observations from the Mi- 
nor Planet Center, The 
ultra-accurate Gaia DR3 observations, combined with the nu- 
merous ground-based observations, represent for the main belt 


 Ihttp: //adams.dm.unipi.it/orbfit/ 
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the equivalent of having very accurate radar measurements for 
the NEOs. Figure[21|shows the accuracy of the along-scan post- 
fit residuals as a projection on the (Aa cos 6, Ad) plane. In the 
projection, two main diagonal lines appear to be more dense: 
they correspond to the initial part of Gaia operations, when an 
ecliptic pole scanning law (EPSL) drove the observations. 


Projection of AAL on Ad (mas) 
Oo 
[o) 
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Fig. 21: Density plot of the projection of the AL post-fit resid- 
uals in the (Aacos 6, Ad) plane. The projection represents the 
quality of Gaia DR3 observations compared to the typical sky- 
plane residuals for ground-based observations. The diagonal 
most prominent lines correspond to the EPSL period. 


Figure also highlights the accuracy of Gaia DR3 data 
when compared to the ground-based post-fit residuals, which are 
typically about 500 mas (two orders of magnitude higher than 
Gaia). 

Gaia DR3 contains 447 NEOs. Of these, 432 can be con- 
sidered as small objects with a diameter smaller than 5 km, and 
197 are small and have a very good orbit uncertainty (the un- 
certainty on the semi-major axis is smaller than 3 x 107° au). 
In the latter list, 24 objects already have a good measurement 
of the Yarkovsky effect in literature. Most of the times, this was 
obtained through radar data. 

We take as an example the case of (3200) Phaethon, the par- 
ent body of the Geminid meteorite shower. Phaethon has al- 
ready a well-established measurement of the Yarkovsky effect. 
From the JPL SBDB A> = —5.56 + 0.68 107! au/d’, where the 
value was obtained using 5090 optical observations and eight 
radar observations. We wish to show that by removing radar ob- 
servations and adding Gaia observations, we are able to find 
a very similar result. The goal is not to prove that we can ne- 
glect radar observations when working with NEOs, but to show 
the strength of Gaia observations when radar is not available, 
for example in the main belt. We used a total of 6723 obser- 
vations. This total includes 356 observations from Gaia DR3. 
We fit the observations to determine the six orbital parameters 
as well as the Az parameter defining the Yarkovsky effect, as 
in|Farnocchia et al,|(2013) and|Del Vigna et al,|2018). We used 
the INPOP10e ephemerides (Fienga et al.!2016) to be consistent 
with the Gaia framework, a gravitational model including the 
eight planets, 23 massive asteroids, and a relativistic model as 
already described in Del Vigna et al.] (2018). We obtain a value 
of Ay = —6.10 + 0.75 au/d*, which is in the lo interval with 
respect to the JPL solution. This result is extremely important 
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because we were able to fit Gaia DR3 with ground-based obser- 
vations, we did not use radar observations, and we did not have 
to manually modify the existing astrometry to obtain a meaning- 
ful result. 

A second example we wish to present is a case for which the 
Yarkovsky effect could not be detected without Gaia astrome- 
try. This is asteroid (1620) Geographos. Using all the available 
observations (5242 optical observations, 105 observations from 
Gaia DR3, and seven radar observations) and the same method 
as described above, we find a value of Ay = —3.25 + 1.01 au/d?. 
It is clear that we are just above the S/N level of 3, and Gaia 
observations allowed the detection. 


6. Photometric performance 


The brightness of an SSO measured at any given epoch depends 
on the observing circumstances. In addition to the distance from 
the Sun and the observer, which can be easily taken into account 
if the orbit of the object is known, they include the rotation of the 
body around its spin axis, and the so-called aspect angle, namely 
the angle between the line of sight of the observer and the ori- 
entation on the celestial sphere of the object spin axis (the as- 
teroid pole). Moreover, the illumination conditions at the epoch 
of observation are critically important. All these parameters vary 
over shorter and longer timescales and mean that any photomet- 
ric measurement of an asteroid is an event that is cannot be re- 
peated in practice. In addition, the measured brightness of an 
SSO observed at any given epoch also depends upon a set of 
constant physical parameters of the object, including its shape, 
spin period, surface geometric albedo, and light-scattering prop- 
erties (Kaasalainen et al./2002). The mechanisms of single and 
multiple scattering of sunlight incident onto the surface of an 
SSO determine the intensity of the flux of scattered sunlight that 
is measured from different directions. These mechanisms depend 
in acomplicated way upon poorly known properties of the object 
surface, including particle sizes, shapes, and optical constants 
(composition), as well as the volume density and surface rough- 
ness of the surface regolith layer. 

For the sake of simplicity, the light-scattering properties of 
an asteroid surface can be described in terms of the dependence 
upon one single parameter, namely the phase angle. This is de- 
fined as the angle between the directions to the Sun and to the 
observer as seen from the target body. 
showed that this can be a reasonable approximation for objects 
with symmetrical shapes and surfaces that scatter incident sun- 
light according to a Lommel-Seeliger surface reflection coeffi- 
cient. It is clear, however, that in the real world, we can expect 
far more complicated situations. 

Based on the above considerations, the validation of sparse 
SSO photometric measurements is not straightforward. Except 
for a handful of objects that were observed in situ by space 
probes, it is practically impossible to make a sufficiently accu- 
rate a priori computation of the expected magnitude of any given 
asteroid observed at an epoch corresponding to given observing 
circumstances. For this reason, the validation procedures devel- 
oped for Gaia DR3 SSO photometric data are based on two basic 
tools: an analysis of the phase-magnitude dependence of SSOs 
in the DR3 catalogue, as explained in Sect. and an analy- 
sis of the test results of inverting sparse photometric data for a 
sample of SSOs, as explained in Sections|6.2]and|6.2.1] The pro- 
cedures and results described in the previous sections represent 
a generalisation and extension of procedures that were adopted 
in the validation of DR2 data (Cellino et al.|2019). In the case of 


DR3, however, the available data are much better in quantitative 
and qualitative respects. 


6.1. Magnitude - phase relation 


Asteroid magnitudes are subject to significant shape-dependent 
periodic variations that are due to rotation around the spin axis. 
In classical ground-based studies, the data normally consist of 
full photometric light curves obtained at different phase an- 
gles. The phase-magnitude relation is then derived by consid- 
ering only magnitude values that were taken at the maximum (or 
mean) brightness value of each light curve in order to avoid to 
mix magnitude data corresponding to different cross sections of 
the rotating body. 

Equally importantly, the photometric data in ground-based 
studies are generally collected during one single apparition of 
an object, namely during a relatively short interval of time (sev- 
eral weeks), during which the object is seen in a nearly con- 
stant geometric configuration. Only the illumination conditions, 
described by the phase angle, vary with time. Recent examples 
of the derivation of phase-magnitude relations are presented by 
and|Mahike etal (2021). 

In the case of Gaia data, however, we deal with sparse mea- 
surements taken during a considerable interval of time and cov- 
ering an interval of phase angles that for main-belt asteroids gen- 
erally ranges from 10° to 30°. The phase-magnitude relation de- 
rived by these data is intrinsically noisy because it includes mea- 
surements taken at epochs corresponding to different illuminated 
cross-sections. This is due both to differences of rotational phase 
around the spin axis and to differences of the orientation of the 
body with respect to the line of sight when data are taken at 
epochs that are sufficiently distant in time. In other words, the 
phase-magnitude curves shown in the Gaia DR3 database are 
contaminated by magnitude variations that are not uniquely due 
to differences in phase angle. This is also a commonly encoun- 
tered situation for the majority of ground-based data, with the 
exception of targeted campaigns. 

These problems can be partially overcome when we consider 
the statistical behaviour of a large number of objects. This is 
made possible by exploiting the fairly large Gaia DR3 database. 
This allows us to adopt some filtering procedures aimed at mim- 
icking the procedures that are traditionally adopted in ground- 
based studies more closely. In particular, of all the available mag- 
nitudes that were measured within the same day, generally corre- 
sponding to two or more consecutive detections in the two FOVs 
of Gaia, only the brightest recorded magnitude was kept in the 
analysis for each object. This was done to limit the noise that is 
uniquely due to the rotation of the object around its spin axis. 
This allowed us to mimic the procedures used in the analysis 
of full light curves taken from the ground more closely, as ex- 
plained above. Moreover, transits for which the apparent Gaia 
magnitude had a nominal error > 0.05 mag were not taken into 
account. It was also decided to discard all objects for which the 
interval of phase angles covered by the observations was exceed- 
ingly narrow, < 5°. As a next step, all asteroids were discarded 
for which the number of accepted magnitude measurements was 
smaller than a limit related to the phase angle interval covered 
by the observations. This limit was set to 4 when the covered 
phase angle interval was larger than 9 deg, it was set to 5 when 
the phase angle interval was between 6 and 9 deg, and it was set 
to 7 when the covered phase angle interval was between 5 and 6 
deg. 

The apparent magnitudes of objects that passed these filters 
were converted into unit distance from the Sun and from Gaia. 
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Fig. 22: Computed phase - Gaia magnitude data for the set of asteroids numbered from 1 to 100 that passed the filtering criteria 
described in the text (top left). The top right panel shows the same, but for objects numbered from1000 to 1100. The bottom 
left panel shows the same, but for asteroids numbered from 10000 to 10100. The bottom right shows the same, but for asteroids 


numbered from 100000 to 100100. 


A linear least-squares fit of magnitude versus phase was finally 
computed. The derived phase - magnitude relations exhibit the 
typical behaviour of asteroids. They are characterised by an over- 
all linear trend in the interval of phase angles covered by Gaia 
observations. Some examples of the obtained phase - magnitude 
relations for different subsets of 100 accepted objects chosen 
in different magnitude regimes (corresponding to objects num- 
bered in different intervals of identification number) are shown 


in Fig. 


The plots shown in Figure correspond to the overlap- 
ping of data corresponding to many different objects, and a gen- 
eral assessment of the statistical behaviour of the whole sam- 
ple of objects passing our selection filter is difficult to derive. 
We therefore performed a more extensive statistical analysis of 
the behaviour of two much larger samples of objects, consid- 
ering their number identifier as a proxy for brightness: the first 
sample included asteroids numbered from 1 to 10000. Of these, 
9107 passed our filtering procedure and were accepted for com- 
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putation of the phase-magnitude relation. The second sample 
was chosen to be that of asteroids numbered from 100,000 to 
110,000. Of these, 4802 were accepted for the phase - magni- 
tude analysis. The decreasing number of accepted asteroids in 
a fainter magnitude regime is related to the decreasing number 
of acceptable transits for increasing faintness. The two samples 
allowed us to compare the phase-magnitude behaviour of ob- 
jects in two clearly different magnitude regimes. Some results 
are shown in Figure[23| They display some general properties of 
the two considered samples. 


We found that the correlation of the obtained linear fits was 
quite variable, with a sharp maximum for values between 0.8 
and 0.9, but far lower values were also included. This was an 
expected consequence of dealing with data that, even after the 
filtering procedures described above cannot completely remove 
the effects caused by the sparseness in time. The nominal er- 
ror of the computed linear slopes was predominantly about 0.01 
mag/degree, but higher values (but very rarely higher than 0.03 
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Fig. 23: Histogram of the computed slopes of the linear fits of the phase - magnitude for asteroids numbered from 1 to 10000 (top 
right). The top right panel shows the same, but for objects numbered from 100000 to 110000. The bottom left panel shows the 
slope error vs slope for asteroids numbered from 1 to 10000 after those with values of the resulting linear correlation < 0.8 and/or a 
number of accepted transits < 8 were removed. The bottom right panel shows the same, but for asteroids numbered from 100000 to 


110000. 


mag/degree) were found in a non-negligible number of cases. A 
histogram of the slopes of the obtained linear fits is shown for the 
two considered samples in the top panels of Figure[23] The vari- 
ety of obtained slopes generally agrees reasonably well with typ- 
ical values mentioned in the literature using ground-based data. 
These typical values range mostly from 0.01 ti 0.04 mag/deg, as 


discussed for instance in|Carbognani et al.|(2019) and/Muinonen| 
(2010). However, some cases of negative slopes (corre- 


sponding to bodies whose brightness would increase for increas- 
ing phase angle, a clearly aberrant result) are found for a mi- 
nority of cases. Some very high values for the slope (> 0.06 
mag/deg) are also found. The existence of these aberrant cases 
is very likely to be due to insufficient removal of transits corre- 
sponding to an exceedingly high variety of observational circum- 
stances, and/or to objects with an insufficient number of avail- 
able measurements. As shown in the bottom panels of Figure[23| 
which show a slope - error (slope) plot for each of the two con- 


sidered samples, but after objects for which the linear correlation 
of the phase - magnitude data turned out to be < 0.08 mag/deg 
and/or the number of accepted magnitude measurements was < 8 
were removed, the resulting linear fits are well compatible with 
typical ground-based values. All the negative values of the slope 
disappear, while the number of high positive values reduces to 
just a few per sample. 


To summarise, we conclude that a preliminary analysis of 
the phase-magnitude relation using the SSO photometric data 
available in Gaia DR3 indicates that in the vast majority of 
cases, the observed phase-magnitude behaviour of the SSOs is 
nicely compatible with the expectations when the unavoidable 
noise arising from the use of limited numbers of sparse photo- 
metric data that are taken in a range of epochs that can corre- 
spond to non-negligible differences in observing circumstances 
are taken into account. We also emphasise that in our analysis, 
we found evidence of differences in the behaviour of objects be- 
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longing to different taxonomic classes, in agreement with cur- 
rent ground-based evidence. In particular, according to[Belskayal 
|& Shevchenko}(2000), dark asteroids tend to have steeper phase- 
magnitude slopes than moderate-albedo asteroids 
[202 1). In this respect, we note that Figure [23] suggests a small 
shift of the peak (from 0.03 to 0.04 mag/deg) in the histogram 
for the objects numbered above 100,000. The reason might be 
that fainter objects tend to be located preferentially in the outer 
region of the asteroid main belt, where dark bodies are predom- 
inant. The distribution of photometric slopes among objects be- 
longing to different taxonomic classes and orbiting at different 
heliocentric distances will be better analysed in future data re- 
leases. 

We are aware that the analysis of the phase-magnitude rela- 
tion described in this section must be considered very prelim- 
inary. Recently, a deep analysis of phase-magnitude data pub- 
lished in Gaia DR2 has been carried out by {Martikainen et al.| 
(2021), who used realistic asteroid shape models and computed 
very good fits of phase-magnitude data against the (H, G,, G2) 
phase function developed by|Muinonen et al.|(2010). We expect 
that a similar analysis will be performed as soon as the Gaia 
DR3 data will become public. 

In this respect, we expect that the future availability of much 
larger numbers of transits per object will allow us in future data 
releases to produce cleaner phase-magnitude plots and to analyse 
the correlations of the photometric behaviour of asteroids with 
their spectroscopic properties and the albedo in more detail and 
to develop a new taxonomy based on Gaia spectroscopic data. 


6.2. Photometric inversion 


Using magnitudes reduced to unit distance, we can analyse the 
variation in the brightness of any given object that is measured at 
different epochs. For a set of sparse photometric measurements 
of the same object, it is convenient to work in terms of bright- 
ness differences with respect to one of the available measure- 
ments (usually the first measurement). In this way, any depen- 
dence upon constant light-scattering properties of the surface can 
be removed, assuming that the surface has homogeneous proper- 
ties. This reduces the time-dependence of the measured bright- 
ness data to a function of the following physical properties: 


1. The rotation period P. This is the rotation of the object 
around its spin axis that continuously modifies the cross sec- 
tion of the illuminated area seen by the observer, depending 
upon the object shape. 


2. The overall shape, which is described by a number of un- 
known parameters. 


3. The orientation of the spin axis with respect to the line of 
sight of the observer (the pole of the object). This corre- 
sponds to two unknown parameters, namely the ecliptic lon- 
gitude and latitude of the pole itself. 


4. The dependence of the brightness upon the illumination con- 
ditions, which is described in terms of the phase angle (see 
(6.1). A simple dependence upon the phase angle summarises 
for the sake of simplicity the overall effect of the mechanisms 
of single and multiple scattering of the sunlight incident onto 
the surface. These mechanisms determine the intensity of the 
flux measured by the observer in different observing circum- 
stances, with a complicated dependence upon poorly known 
surface properties, including albedo, texture, and roughness. 
These can be assumed to be a constant but unknown function 
of the phase angle for any given object. 
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Based on these considerations, it is in principle possible to 
develop numerical codes to determine the unknown physical pa- 
rameters of an object whose brightness has been measured by 
Gaia in a sufficiently large number of observed transits. A code 
like this, based on a genetic algorithm, has been developed for 
the purposes of Gaia data processing. It will be used to produce 
results of SSO photometric inversions in future Gaia data re- 
leases. The algorithm assumes that the shapes of the object are 
triaxial ellipsoids, described by two parameters (axial ratios), 
and that there is a linear variation in magnitude as a function 
of the phase angle (see (6. 1p. This algorithm has already been 
adopted in preliminary analyses of Gaia DR2 data 
(2019), and has been used to validate Gaia DR3 photometric data, 
as explained in what follows. 

A preliminary step was identifying objects with reliable pre- 
dictions of the brightness at any given epoch of observation. 
We profit in this way from a detailed knowledge of the phys- 
ical properties characterising the object. In principle, the best 
possible validation test must consist of comparisons between ex- 
pected and measured magnitudes for a set of objects for which 
our knowledge of their physical parameters is extraordinarily ac- 
curate, being based on in situ measurements carried out by space 
probes. The number of these objects is unfortunately extremely 
small, and they deserve a separate treatment. It is therefore nec- 
essary to take advantage of larger data sets of ground-based as- 
teroid photometric data. Decades of ground-based photometry, 
mostly at visible wavelengths, have produced large catalogues 
of asteroid light curves. The rotation periods of about 10 000 of 
these asteroids have been derived with good accuracy. A smaller 
number of these objects have been observed in a variety of ob- 
servation circumstances sufficient to derive accurate estimates of 
the spin axis direction for them. In many cases, more than one 
pole solution is found to be compatible with available data. In 
a large number of cases, an overall shape, derived using com- 
plex algorithms of light-curve inversion, has also been obtained. 
Shape details, however, are of limited importance in our case, 
because we compute photometric inversion using a simple triax- 
ial ellipsoid model, which is only a first approximation of what 
the real shape of an asteroid can be. In our analysis, we used cur- 
rently available catalogues of asteroid rotation periods and pole 
coordinates. In particular, we massively exploited the database 
of asteroid models from i version paid us (DAMIT), which 
is publicly available (Durech et al.| (Durech et al.| 2010)]" | For objects not in- 
cluded in DAMIT, we took the sae from the asteroid light 
curve database catalog, available at the web site of the NASA 
Planetary Data System. 


We limited our analysis to the brightest asteroids, numbered 
from 1 to 1000, for which at least 30 accepted transits were avail- 
able in the Gaia DR3 database. In particular, we accepted only 
transits for which the nominal error of the G magnitude was not 
larger than 0.02 mag because the results of photometric inversion 
can be negatively affected when data of insufficient quality are 
used. In any case, our sample of asteroids includes objects that 
are sufficiently bright to have systematically smaller magnitude 
uncertainties than the above limit. The resulting sample used for 
our photometric inversion test includes 430 asteroids. 


We note that the results of our photometric inversion algo- 
rithm are also able to distinguish between prograde and retro- 
grade rotation, according to IAU criteria. In particular, objects 
found to have a retrograde rotation have negative assigned val- 
ues of their estimated P value, while the determined spin axis di- 


T \https://astro.troja.mff.cuni.cz/projects/damit/ 
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Fig. 24: Relative difference between the absolute value of the 
difference between the resulting P solution and the P value de- 
termined by ground-based observations in units of the ground- 
based P value for each object of the considered sample. 


rection, expressed in ecliptic longitude and latitude, has always 
an assigned positive value for the latitude of the pole. 


It is important to note that the adopted algorithm for photo- 
metric inversion produces for each object a set of 15 different 
inversion solutions for each object. The reason is that an intrin- 
sic property of the adopted genetic algorithm is that it does not 
always converge to a unique final solution. This happens because 
the evolution of the population of possible solution parameters 
can often evolve along a dead branch, leading to a poor-quality 
solution. Fifteen different genetic solutions have been found to 
be a reasonable compromise between the need to obtain a good- 
quality solution and the need of minimising the CPU execution 
time. 


We stress that at this stage, only the best obtained solution 
for each object, namely the solution that produces the smallest 
residuals with respect to the measured magnitudes at the avail- 
able transits, was considered and compared with ground-based 
results. This choice has some consequences, because in several 
cases, more than one inversion solution gave equivalent resid- 
uals. We defined as equivalent residuals those that differ by no 
more than 0.0015 mag. As a consequence, the inversion solution 
is not unique in some cases. In these cases, the nominally best 
solution may not correspond to the rotation period and pole listed 
in ground-based catalogues, whereas another equivalent solution 
that is not considered for the moment would correspond to these 
values. For this reason, the results presented here are conserva- 
tive. 

We recall that for the purposes of inversion of a set of sparse 
photometric measurements, it is of paramount importance to 
have a good sampling of the possible observing circumstances 
for any given object. This means that the data should include 
measurements that adequately sample the whole interval of 360 
deg in ecliptic longitude. This is especially important for the de- 


termination of the spin axis orientation (the asteroid pole). How- 
ever, this is not yet the case for the data that are available in 
Gaia DR3. For each object, large gaps exist in the interval of 
the covered ecliptic longitudes. The publication of photometric 
inversion of Gaia data is scheduled as an end-of-mission task. 
The results of the preliminary inversion attempts presented here 
must be considered as no more than a useful tool for the sci- 
entific validation of Gaia DR3 data. We expect to obtain much 
better inversion solutions in future data releases as new measure- 
ments will become available. 

We considered as a successful determination of the rotation 
period P an inversion solution for which the absolute value of the 
difference between the resulting P solution and the P value de- 
termined by ground-based observations, expressed in units of the 
absolute value of the ground-based P value, expressed in hours, 
is not higher than 0.001. This criterion takes into account the fact 
that for fast rotations of just a few hours, small errors in P lead to 
strong differences in the rotational phase of the object, which is 
computed at epochs differing by few years. On the other hand, in 
the case of very long rotation periods, longer than several dozen 
or hundreds of hours, it is not reasonable to impose a required 
accuracy of about a few seconds, for example, on the determina- 
tion of the period. 

Based on our adopted criterion, we obtained the correct P 
solution for 229 out of 430 asteroids of our sample. These re- 
sults are shown in Fig. [24] This figure shows an interesting fea- 
ture: in addition to the 229 cases of correct P determinations, 
a significant number of cases exist,for which the P value deter- 
mined by photometric inversion is nearly exactly twice the P 
value determined from ground-based photometry. This feature 
is not entirely unexpected considering the assumptions of the 
adopted model. In particular, the algorithm assumes that over a 
full rotation of the object, the brightness reaches two maxima 
and two minima, which are expected to be (nearly) equal. In 
the real world, however, and in particular when an object is not 
strongly elongated and the maxima and minima tend to be shal- 
low and/or strongly asymmetric, the inversion algorithm might 
derive a rotation period that can be twice the correct one. More- 
over, at least in some cases, a derived double period might be 
indicative of a photometric behaviour that is dominated not by 
the shape, as assumed by the inversion model, but by variation 
in surface albedo. This is the case, for instance, of the large aster- 
oid (4) Vesta (Cellino et al.]1987), which is known to have a light 
curve producing only one maximum and minimum per cycle. 

It seems therefore that the results shown in Fig. [24] are very 
encouraging when the limits of the adopted shape model, the 
small minimum number of accepted transits per object, the still 
limited variety of observing circumstances, and the conservative 
criteria of definition of the inversion solution are taken into ac- 
count. Ambiguous cases are not considered, together with the the 
non-negligible number of solutions corresponding to a P twice 
as large as ground-based determination. 

To determine the pole, the situation is intrinsically more 
complicated. In the vast majority of cases, two or more different 
pole solutions per object are listed in the literature. Moreover, in 
our assumption of a triaxial ellipsoid shape model, an ambiguity 
of 180 degrees in the determination of the ecliptic longitude of 
the pole may be present in some cases due to the symmetry of 
this shape model. In particular, a 180° ambiguity on the longi- 
tude of the pole can be triggered by a low orbital inclination of 
the object or by an unfavourable distribution of the observations 
in ecliptic longitude. Moreover, simulations have shown that a 
low ecliptic latitude of the pole of the object tends to make the 
photometric inversion more challenging with the assumed shape 
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Fig. 25: Comparison of Gaia DR3 inversion solutions for the spin period P and ground-based P determinations for all asteroids of 
our sample with a pole quality code = 2 (best inversion solutions; top panel). The top right panel shows the same, but for objects 
with a pole quality code =1. The bottom panels show the same as the top panels, but for objects with a quality code solution = 0 
(corresponding to an unknown pole solution from ground-based data) and for objects with a pole quality code = -1 (corresponding 
to complete disagreement with any existing ground-based pole solution). The upper line represents periods derived from Gaia DR3 


photometry that are exactly twice those obtained from the ground. 


model (Santana-Ros et al.|/2015). We limited our analysis to a 


comparison between the obtained pole solutions and the DAMIT 
database of asteroid poles, which is considered to be the most 
accurate list of ground-based asteroid period and pole determi- 
nations. Based on the results of our analysis, we decided to sub- 
jectively define a small number of pole solution quality classes 
(QC). In particular, we assigned pole QC = 2 to pole solutions 
that differed by no more than about 10 deg (separately) in eclip- 
tic longitude and latitude with respect to one existing DAMIT 
pole solution, and we also imposed that the rotation period de- 
termination was correct. We assigned pole QC = 1 to objects for 
which the best pole solution was not so close to a DAMIT pole 
solution, independently of the obtained period. We assigned pole 
QC = -1 in situations in which the obtained pole solution had lit- 
tle to do with any existing DAMIT solution. Finally, we assigned 
pole QC = 0 to objects for which no DAMIT pole solution exists. 
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Figure[25]shows plots of the period determined by inversion 
of Gaia DR3 data versus the P value known from ground-based 
determinations for objects belonging to different pole solution 
QCs. The figure shows that except for the case of pole solution 
QC = 2, for which a (nearly) perfect agreement with ground- 
based period solutions is imposed by definition, when objects 
belong to different pole solution QCs, the relation between the 
QC of the pole solution and the success in the period determina- 
tion by inversion of Gaia DR3 data is not always obvious. In par- 
ticular, poor pole determinations (red points in the bottom right 
panel of Fig. do not correspond to a larger fraction of erro- 
neous period determinations. We note, however, that the fraction 
of cases in which the period from Gaia DR3 photometric data 
inversion tends to disagree with ground-based period determina- 
tions tends to increase when we consider objects for which no 
pole solution from ground-based data is available (see the bot- 
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Table 4: Rotation period (P, in hours), the ecliptic longitude and 
latitude of the pole (A and £, in degrees), and the Tholen taxo- 
nomical classes for asteroids (21) Lutetia and (2867) Steins. For 
the latter, there are two possible pole solutions I and Il. Num- 
bers in parentheses depict the published uncertainty in units of 
the last digit shown. 


(2867) Steins 
6.0468 1(2) 


(21) Lutetia 
P 8.168270(1) 


A,B 52.2(4), -7.8(4) I: 96(5), -85(5) 
II: 142(5), -83(5) 
Class M E 


tom left panel in Fig.[25). This might be considered as an indica- 
tion that asteroids for which no ground-based pole determination 
is available may well be challenging cases, for which even the 
ground-based period solutions could be more uncertain and pos- 
sibly incorrect. The relatively high number of objects with long 
rotation periods among those determined by ground-based data 
should be noted in particular, and for which ground-based data 
are so far insufficient to derive a pole solution. Long rotation pe- 
riods correspond in many cases to relatively large uncertainties 
in the determination of the period. 

Based on the results shown in this and the previous section, 
we conclude that Gaia DR3 photometric data are of a good if 
not excellent quality because we successfully produced a correct 
inversion of a large number of objects using a clearly simplis- 
tic shape model and relatively small numbers of measurements 
covering a still partial fraction of the possible observing circum- 
stances. We limited our analysis for the moment to measured 
Gaia DR3 magnitudes with nominal errors not exceeding 0.02 
mag. The results of our photometric inversion attempts seem to 
be very encouraging. In principle, we cannot rule out the possi- 
bility that in some cases, the error bar of some Gaia DR3 magni- 
tude might be higher than the nominal value. In some cases, this 
might affect the results of photometric inversion negatively, and 
might explain some incorrect results. 

More convincing tests should be based on the analysis of 
data of asteroids for which our knowledge of the rotation period 
and spin axis direction are of the best possible reliability. This is 
the case of a very small number of objects that were visited in 
situ by space probes. The results for some of them, for which we 
have a reasonable number of Gaia DR3 observations, are shown 
in the next section. 


6.2.1. Photometry of (21) Lutetia and (2867) Steins 


Our validation consists of comparing the observed Gaia pho- 
tometry to photometry computed for known asteroids that have 
accurately determined rotation periods, pole orientations, and 
high-resolution shape models. Gaia DR3 contains photometric 
data for more than a dozen asteroids studied by space missions. 
By using the shape models, rotational parameters, and taxonom- 
ical classifications of (21) Lutetia and (2867) Steins, which are 
asteroids that were visited by the ESA Rosetta space mission, 
we studied whether it is possible to validate the Gaia DR3 SSO 
photometry. We note that (21) Lutetia and (2867) Steins were 
assessed earlier in the documentation of Gaia DR2. 

For (21) Lutetia and (2867) Steins, the Planetary Data Sys- 
tem provides the shape models illustrated in Fig.[26]as well as the 


rotation periods and pole orientations described in Table|4|(Farn- 
2011). The table also includes the Tholen taxonomical classes of 
the two asteroids. 


The Gaia DR3 photometric measurements for these two as- 
teroids are plotted as a function of the observation epochs, ex- 
pressed in days after the first observation, and against the phase 
angle in degrees, in Fig. Both phase-magnitude relations 
show a decreasing trend in disk-integrated brightness with in- 
creasing phase angle. Furthermore, the apparent slope of de- 
creasing brightness is steeper for Lutetia, in agreement with 
Lutetia and Steins being lower-albedo M-type and higher-albedo 
E-type objects in the Tholen taxonomy, respectively. 


We studied (21) Lutetia with the 23 transit magnitudes ob- 
tained by Gaia (case I) using a shape model derived from com- 
bined ground-based relative photometry of 50 light curves with 
4012 photometric points in total (Durech et al.!2010) and Gaia 
photometry with convex inversion methods (Muinonen et al. 
The results showed that the Gaia measurements contributed to 
the shape modelling with reasonable residuals, with a low RMS 
value of 0.016 mag. The RMS value for the entire data set 
of 51 light curves was 0.015 mag. Next (case II), it was seen 
that using the best available shape model (Farnham]2013), con- 
structed by combining a high-resolution shape model based on 
disk-resolved imaging by Rosetta and a lower-resolution model 
based on ground-based relative photometry and silhouette ob- 
servations with adaptive optics, it was not possible to produce 
a straightforward fit to the Gaia observations (see Fig. |28). By 
optimising the rotational phase of Lutetia and the slope of the 
phase curve implied by the Lommel-Seeliger scattering model, 
the high-resolution shape model resampled at a 5-degree resolu- 
tion reproduced the Gaia photometry with a high RMS value of 
0.047 mag, whereas for the entire data set, the RMS value was 
0.017 mag. We interpreted the RMS values for cases I and II as 
giving a strong indication that the Gaia data offer new informa- 
tion on the shape or surface properties of Lutetia. 


A closer inspection of the time distribution of the 23 Gaia 
DR3 transit magnitudes shows that for 12 of them, the Lutetia 
hemisphere that predominantly visible to Gaia was the hemi- 
sphere that was not mapped by Rosetta. It is therefore highly 
probable that the Gaia data provide additional information about 
this portion of the asteroid surface. In particular, the Gaia 
photometry, with its absolute phase angle dependence, relates 
the size and geometric albedo characteristics of the two hemi- 
spheres. 


We used the high-resolution shape resampled in 5-degree 
resolution to model the photometry of the remaining 11 observa- 
tions that were obtained with the hemisphere mapped by Rosetta 
in view directly (case III). Allowing for optimisation in only 
the rotational phase and the slope of the phase curve implied 
by the Lommel-Seeliger scattering model, the high-resolution 
shape model reproduced the Gaia photometry with an excellent 
RMS value of 0.019 mag (Fig. [28). We then applied an analo- 
gous approach to the 12 observations obtained with the hemi- 
sphere that was not mapped by Rosetta (case IV). With con- 
stant geometric albedo characteristics, the RMS value was high 
at 0.057 mag. 


The photometric slopes implied by the Lommel-Seeliger 
scattering model were then compared in detail. In the cases 
of convex optimisation using all Gaia measurements (case I), 
optimisation of the rotational phase and slope using the 
high-resolution shape model and all Gaia measurements (II), 
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Fig. 26: Shape model of asteroid (21) Lutetia (left). Shape model of asteroid (2867) Steins (right). The models are based on in situ 
images obtained by the Rosetta space mission. In the panels, the top and bottom plates on the left correspond to polar views along 
the z-axis (axis of rotation). The top (bottom) plates in the middle and to the right correspond to viewing along the x-axis (y-axis). 
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Fig. 27: Observed G magnitudes of asteroid (21) Lutetia included in Gaia DR3, relative to the mean of the magnitudes (Go), as a 
function of days after the first observation (fo) (top left). The top right panel shows the same as the left panel for (2867) Steins. The 
bottom panel shows the same as the top panel, but the magnitudes are depicted against the phase angle. 


measurements corresponding to the hemisphere observed by 
Rosetta (III), and measurements corresponding to the hemi- 
sphere that was not observed by Rosetta (IV), the slopes 
were 1.89 mag rad7!, 1.56 mag rad~!, 1.65 mag rad7!, and 
1.47 mag rad7!, respectively. First, case I stands out as show- 
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ing the steepest photometric slope for a convex shape solution. 
A plausible explanation is that global self-shadowing effects 
due to a non-convex shape steepen the photometric slope. Sec- 
ond, case IV stands out as showing the shallowest photomet- 
ric slope for the high-resolution non-convex shape model. Third, 
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Fig. 28: Observed G magnitudes of asteroid (21) Lutetia included in Gaia DR3 as a function of days after the first observation (fo, 
top left) and against the phase angle (a, top right) together with the modelled magnitudes. The bottom left panel shows the same 
as in top right panel, but only for the observations concerning the Lutetia hemisphere observed by Rosetta. The bottom right panel 
shows the same as in the top right panel, but only for the observations of the Lutetia hemisphere that was not observed by Rosetta. 
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Fig. 29: Observed G magnitudes of asteroid (2867) Steins included in Gaia DR3 against the phase angle a together with the 
modelled magnitudes. Two models are shown, corresponding to the pole orientations I (left) and I (right) given in Table|4] 


it is beyond the scope of the present study to improve the Lutetia 
model, the present analysis indicates varying scattering proper- 
ties for the hemisphere that was not imaged by Rosetta. This is 
also supported by the fact that all the ground-based light curves, 


case III shows the most realistic photometric slope for the high- 
resolution shape model: the least-squares fit to the Gaia data is 
successful with a small RMS-value. Finally, the case II photo- 
metric slope is a compromise between cases III and IV. Whereas 
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consisting of extensive observations of both hemispheres and 
treated in relative sense in the validation, were well fitted by the 
high-resolution shape model. 

The final photometric slope analysis for Lutetia was carried 
out in case III using Markov chain Monte Carlo (MCMC) sam- 
pling. The photometric slope for the hemisphere observed by 
Rosetta was estimated to be (1.67740.075) mag rad7! using the 
high-resolution shape model resampled at a 5-degree resolution. 
Based on[Penttilé et al,|(2016) [Muinonen etal, (2020), and|Mar- 
(2021), these photometric slopes agree with those 
for Tholen S- and M-class asteroids. 

For Lutetia, the apparent G magnitudes range from 11.1 to 
13.1 mag, with minute nominal observational errors of 0.00033 — 
0.00070 mag. Our RMS values above are significantly higher. 
We conclude that Gaia photometry has an accuracy better than 
0.01-0.02 mag with the limitations on first, the shape model ac- 
curacy, which does not allow us to push the analysis further. Sec- 
ond, we assume that the same Lommel-Seeliger scattering model 
is valid everywhere on the surface of Lutetia. In particular, we 
assume that a single combination of geometric albedo and scat- 
tering phase function is representative of the entire surface of 
Lutetia. Third, it is not possible to exclude that some anomalous 
brightness values are present. 


The results obtained for (2867) Steins, whose high- 
resolution shape model is not limited by hemisphere scales sim- 
ilar to those of (21) Lutetia, support the idea that Gaia photom- 
etry is indeed accurate. In the case of this object, there are two 
pole solutions (Steins cases I-II, see Table|4) that essentially only 
differ in terms of the zero-point of the rotational phase. By di- 
rectly using the shape model to reproduce the Gaia data, resam- 
pled at a 5-degree resolution, the RMS values of the observed- 
computed magnitudes are 0.022 mag and 0.023 mag for the two 
pole solutions. These fits are based on 30 photometric points, 
3 of which had to be omitted as outliers. Again, the fits were 
obtained by optimising only the rotational phase and the photo- 
metric slope implied by the Lommel-Seeliger scattering model. 


We continued the analysis of Steins photometry by resam- 
pling the shape model at a higher 3-degree resolution and si- 
multaneously fitting five parameters, that is, the rotation period, 
pole orientation, rotational phase, and photometric slope. The 
RMS values for the two pole solutions were lower at 0.018 mag 
and 0.019 mag, respectively. To compute these fits, the rotational 
parameters were regularised to lie within their domains of un- 
certainty (Table (4). The photometric slope obtained values of 
1.515 mag rad~! (case I) and 1.522 mag rad~! (II), that is, the 
values were essentially equal. 

Finally, MCMC sampling was carried out for the rota- 
tional phase and photometric slope in cases I-II using the high- 
resolution shape model re-sampled at 5-degree resolution. For 
5000 samples, the means and standard deviations of the pho- 
tometric slope obtained the values of (1.497+0.064) mag rad7! 
(case I) and (1.458+0.062) mag rad7! (II). The values agree mu- 
tually, and they are realistic for high-albedo E-class asteroids 
(Pentil et al.|2016) 
2021). The mean values agreed within the given uncertainties 
with the least-squares values for the resampled 5-degree and 3- 
degree shape models. 


For Steins, the apparent G magnitudes ranged over 16.3 — 
18.4 mag, with nominal observational errors within 0.0028 — 
0.014 mag. Considering the simplicity of the surface-scattering 
modelling we used, our results can be considered completely 
satisfactory (see Fig. [29). The remaining limitations in the case 
of (2867) Steins are related to details of the shape model and 
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surface-scattering characteristics. In particular, certain assump- 
tions were made on the scattering properties when the high- 
resolution shape model was derived from the Rosetta images. 
It remains possible that the three observations that we omitted 
as outliers were omitted due to modelling issues rather than low 
observational accuracy. 

In conclusion, we cannot rule out the possibility that the sam- 
ple published in Gaia DR3 could still include a non-negligible 
fraction of anomalous data. We recommend detailed analyses 
and careful verifications when the Gaia DR3 photometry of as- 
teroids is applied. However, especially for objects as bright as 
Lutetia, our current approach to magnitude prediction cannot in- 
vestigate Gaia photometry of SSOs at the level that would be 
required by the very high accuracy. 


7. Conclusions 


The Solar System data processing, developed and trained over 
several years, reaches maturity with Gaia DR3 by providing the 
expected large survey for asteroids and planetary satellites. The 
different data types (astrometry, photometry, and low-resolution 
reflectance spectra) and their accuracy and homogeneity are an 
impressive achievement of the Gaia mission. We have summa- 
rized the approach that we followed for the processing pipeline 
and in particular explained the procedures that clean the data set 
and provide a reliable outcome. 

The highly accurate epoch astrometry obtained by Gaia is 
confirmed by the orbit adjustment, resulting in residuals at sub- 
milliarcsecond level for G < 18. The astrometric performance is 
clearly improved with respect to Gaia DR2. 

This accuracy for a large number of asteroids brings new ca- 
pabilities of investigation, revealing effects related to the par- 
tially resolved shape of asteroids. The amplitude of the photo- 
centre wobble that we measured on (21) Lutetia shows that this 
effect probably affects the orbit computation even for asteroids 
that are three to four times smaller. The capability of detecting 
the wobbling, associated with the orbital motion of a satellite, 
discloses an impressive domain of investigation for the search 
of asteroid binaries with the astrometric method. This approach 
was simply not possible before. A comprehensive exploration of 
Gaia DR3 should reveal a variety of binary systems. 

A search for the best orbital modelling for near-Earth objects 
that can exhibit the Yarkovsky drift shows that in some circum- 
stances, Gaia can detect this effect, even in the absence of the 
radar ranging data that have been essential in the recent past. 
A careful coupling to ground-based data will fully disclose this 
potential. 

The G-band photometric data are the other valuable source 
of information that will be exploited to obtain new constraints on 
the rotation and the shape of a very large number of objects. For 
them, it remains difficult to assess the quality at the level of their 
expected uncertainty, simply because there is no direct compar- 
ison to other data sets of comparable accuracy. The simulated 
photometry from accurate shape models brings, however, pos- 
itive results, despite the remaining uncertainties. It also shows 
that some outliers remain probably present among the released 
measurements. 

In conclusion, our review shows that the data processing of 
Solar System data in Gaia DR3 has produced an extremely rich 
data set standing out for its unique properties with respect to 
other existing surveys in many aspects. Gaia DR3 will certainly 
be exploited in many ways by the community of planetary sci- 
entists. New properties and features beyond those illustrated in 
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this article will probably be found, and will be a major driver for 
improvements of the data quality in future data releases. 
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